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Ph' Abstract 

P5 I We revise fundamental concepts in the dynamics of open quantum 

^ ■ systems in the hght of modern developments in the field. Our aim is to 

qh. present a unified approach to the quantum evolution of open systems 

that incorporates the concepts and methods traditionally employed by 
' different communities. We present in some detail the mathematical 

■ structure and the general properties of the dynamical maps underlying 
open system dynamics. We also discuss the microscopic derivation of 

(SJ ' dynamical equations, including both Markovian and non-Markovian 

■ evolutions. 
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1 Introduction 

To write an introduction to the dynamics of open quantum systems may seem 
at first a complicated, albeit perhaps unnecessary, task. On the one hand, the 
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field is quite broad and encompasses many different topics which are covered 
by several books and reviews On the other hand, the approach taken 

to study the dynamics of an open quantum system can be quite different 
depending of the specific system being analyzed; for instance, references in 
quantum optics (e.g. [T|[T0 | [TT | [T5 | [T7 t [26] ) use tools and approximations which 
are usually quite different from the studies focused on condensed matter or 
open systems of relevance in chemical physics (e.g. [HI [131 EH EDI E3])- More- 
over, very rigorous studies by the mathematical physics community were car- 
ried out in the 1970's, usually from a statistical physics point of view. More 
recently, and mainly motivated by the development of quantum information 
science, there has been a strong revival in the study of open many body sys- 
tems aimed at further understanding the impact of decoherence phenomena 
on quantum information protocols |38j . 

This plethora of results and approaches can be quite confusing to novices 
in the field. Even within the more experienced scientists, discussions on 
fundamental properties (e.g. the concept of Markovianity) have often proven 
not to be straightforward when researches from different communities are 
involved. 

Our motivation in writing this work has been to try to put black on white 
the results of several thoughts and discussions on this topic with scientists of 
different backgrounds during the last three years, as well as putting several 
concepts in the context of modern developments in the field. Given the ex- 
tension of the topic and the availability of many excellent reviews covering 
different aspects of the physics of open quantum systems, we will focus here 
on one specific issue and that is embedding the theory as usually explained 
in the quantum optical literature within the mathematical core developed in 
the mathematical physics community. Finally, we make a connection with 
the current view point of the dynamics of open quantum systems in the light 
of quantum information science, given a different, complementary, perspec- 
tive to the meaning of fundamental concepts such that complete positivity, 
dynamical semigroups, Markovianity, master equations, etc. 

In order to do this, we will restrict the proof of important mathemat- 
ical results to the case of finite dimensional systems; otherwise this study 
would be much more extensive and potentially interested people who are not 
expert in functional analysis or operators algebras, could not easily bene- 
fit from it. However, since the most important application of the theory of 
open quantum system nowadays is arguably focused on the field of controlled 
quantum technologies, including quantum computation, quantum communi- 
cation, quantum metrology, etc., where the most useful systems are finite (or 
special cases of infinite systems), the formal loss of generality is justified in 
this case. The necessary background for this work is therefore reduced to 
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familiarity with the basics of quantum mechanics and quantum information 
theory. 

Furthermore, we have tried to incorporate some novel concepts and tech- 
niques developed in the last few years to address old problems such that those 
exemplified by the use of dynamical maps, Markovian and non-Markovian 
evolutions or microscopic derivations of reduced dynamics. We have orga- 
nized the presentation as follows. In section [2] we introduce the mathematical 
tools required in order to ensure a minimum degree of self-consistency in the 
study; a mathematical experienced reader may skip this part. In section |3] we 
review very succinctly the concept of "closed quantum system" and discuss 
its time evolution. We explain in section H] the general features of dynamics 
of open quantum systems and its mathematical properties. Section deals 
with the mathematical structure of a very special important dynamics in 
open quantum systems, which are the Markovian quantum processes. The 
microscopic derivation of Markovian quantum dynamics is explained in sec- 
tion [6l whereas some methods for obtaining non-Markovian dynamics are 
briefly introduced in section [71 Finally, we conclude the review by making a 
summary of the main ideas we have discussed, and the topics that we have 
left out, in the section |8l 



2 Mathematical tools 



2.1 Banach spaces, norms and linear operators 

In different parts within this work, we will make use of some properties of 
Banach spaces. For our interest it is sufficient to restrict our analysis to finite 
dimensional spaces. In what follows we revise a few basic concepts [391 - HTj . 

Definition 2.1. A Banach space ^ is a complete linear space with a norm 

Recall that a space is complete if every Cauchy sequence converges in it. 
An example of Banach space is provided by the set of self-adjoint trace-class 
linear operators on a Hilbert space "H. This is, the self-adjoint operators 
whose trace norm is finite || ■ ||i, 

||A||i = Tr = Tr A G «B. 

Now we focus our attention onto possible linear transformations on a 
Banach space, that we will denote by T(Q3) : ® i— t- The set of all of them 
forms the dual space 
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Proposition 2.1. The dual space of^, 23* is a (finite) Banach space with 
the induced norm 

\\T\\= sup ^ri^= sup ||T(x)||, Tg*B*. 

Proof: The proof is basic. 

Apart from tlie triangle inequality, the induced norm fulfills some other 
useful inequalities. For example, immediately from its definition one obtains 
that 

< ||r||||a;||, Vxe^B. (1) 

Moreover, 

Proposition 2.2. Let Ti and T2 be two linear operators in 03*, then 

IIT1T2II < IIT1IIIIT2II. (2) 

Proof. This is just a consequence of the inequality ([1]); indeed, 

\\T^T2\\= sup ||TiT2(x)|| < sup ||Ti|| ||T2(x)|| = ||Ti|| ||T2||. 
xe'S,||x||=i xe'B,||x||=i 

Q.E.D. 

The next concept will appear quite often in the forthcoming sections. 

Definition 2.2. A linear operator T on a Banach space 03, is said to be a 
contraction if 

\\T{x)\\ < \\x\\, Vx e 03, 

this is ||T|| < 1. 

Finally, we briefly revise some additional basic concepts that will be fea- 
tured in subsequent sections. 

Definition 2.3. For any linear operator T on a (finite) Banach space !B, its 
resolvent operator is defined as 

Rt{X) = {Xl-T)-\Xe C. 

It is said that A G C belongs to the (point) spectrum of T, A G ciT), if -Rt(A) 
does not exist; otherwise A G C belongs to the resolvent set of T, A G g{T). 
As a result (j(T) and g(T) are said to be complementary sets. 

The numbers A G cr(T) are called eigenvalues of T, and each element in 
the kernel x G ker(Al— T), which is obviously different from {0} as (Al— T)~^ 
does not exist, is called the eigenvector associated with the eigenvalue A. 
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2.2 Exponential of an operator 

The operation of exponentiation of an operator is of special importance in 
relation to linear differential equations. 

Definition 2.4. The exponential of a linear operator L on a finite Banach 
space is defined as 

n=0 

The following proposition shows that this definition is indeed meaningful. 
Proposition 2.3. The series defining is absolutely convergent. 
Proof. The following chain of inequalities holds 



n=0 



^ IlLr 11,11 

^ n! ^ n\ 

n=0 n=0 



Since e"^" is a numerical series which is convergent, the norm of is upper 
bounded and therefore the series converges absolutely. 

Q.E.D. 

Proposition 2.4. // Li and L2 commute with each other, then e^^i"'"^^) = 

Proof. By using the binomial formula as in the case of a numerical series, we 
obtain 

' k T n~k 



n=0 n=0 fc=0 ^ ^ 

00 n ^ 00 00 



^^k\in-k)\ ^ ^ p\q\ ^ ^ 

n=0 fc=0 ^ ' p=0 q=0 ^ ^ 

Q.E.D. 

The above result does not hold when Li and L2 are not commuting op- 
erators, in that case there is a formula which is sometimes useful. 

Theorem 2.1 (Lie- Trotter product formula). Let Li and L2 be linear oper- 
ators on a finite Banach space, then: 



n— >oo V 
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Proof. We present a two-step proof inspired by [12]. First let us prove that 
for any two linear operators A and B the following identity holds: 



n-l 



A" - = ^ A''{A - B)B''-''-\ (4) 

fc=0 

We observe that, in the end, we are only adding and subtracting the products 
of powers of A and B in the right hand side of (|4j), so that 

A" - = A" + {A""^B - - B"" 

= A'' + {A^'-^B - A''-^B + A^'-^B^ - A^-'^B'^) - B"" = 
= A" + {A^'-^B - A^'^^B + _ ^ _ _ ^ 

■ ■ ■ + AB"~^ - AB""-^) - B". 

Now we can add terms with positive sign to find: 

n.-l 

and all terms with negative sign to obtain: 

-A^'-^B - ... - - 5" = - ^ A'^'fi"' 



/I £? 

fc=0 

Finally, by adding both contributions one gets that 

n— 1 n— 1 

fc=0 fc=0 

as we wanted to prove. 
Now, we denote 

and take A = Xn and = y„ in formula (jlj), 

n-l 
fc=0 

On the one hand, 

= II (^e(''^i+'^2)/nj'= II <; ||g(Li+L2)/n||fc < g||Li+L2 ||fc/n 

<e(ll^ill+ll^2||KVn<ell^ill+l!i2||^ VA; < n. 
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Similarly, one proves that 
So we get 



n-l 



n-l 



k=0 k=0 
n-l 



k=0 



On the other hand, by expanding the exponential we find that 



lim ?i||X„— = lim n 

n—^oo n—^oo 



^ fci I n / ^ k\j\ \ n J \n 



k=0 



< lim n 

n— >oo 



' k,j=0 



0. 



Therefore we finally conclude that: 



lim 

n— >oo 



lim ||x:-i;"||=o. 



Q.E.D. 



Using a similar procedure one can also show that for a sum of operators, 
the following identity holds: 



N 



lim TT 



e " 



(5) 



2.3 Semigroups of operators 

Definition 2.5 (Semigroup). A family of linear operators Tt {t > 0) on a 
finite Banach Space forms a one-parameter semigroup if 

1. TtTs = Tt+s, Vt, s 

2. To = 1. 
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Definition 2.6 (Uniformly Continuous Semigroup). A one-parameter semi- 
group Tf is said to be uniformly continuous if the map 

t^Tt 

is continuous, this is, hm^.^^ ||Tj — T^H = 0, Vs. This kind of continuity is 
normally referred to as "continuity in the uniform operator topology" |39] . 
and it is sufficient to analyze the finite dimensional case, which is the focus 
of this work. 

Theorem 2.2. IfTt forms a uniformly continuous one-parameter semigroup, 
then the map t ^ Tt is differentiable, and the derivative of Tt is given by 

dTt 

-r- = LTt, 
dt 

with L = ^\t=o- 

Proof. Since Tt is uniformly continuous on t, the function V{t) defined by 

V{t) = [ Tsds, t > 0, 
Jo 

is differentiable with = Tt (to justify this, it is enough to consider 

the well-known arguments used for the case of R, see for example |13] , but 
substituting the absolute value by the concept of norm). In particular. 



t^o t t^o t dt 



= To = 1, 

t=o 



this implies that there exist some to > small enough such that V{tQ) is 
invertible. Then, 



Tt = V~\t,)V{t^)Tt = V-\to) r TsTtds = V-\to) r Tt+sds 

Jo Jo 

rt+h 

= ^~'(to)J^ Tsds = V-\to)[V{t + to)-V{t)], 

since the difference of differentiable functions is differentiable, Tt is differen- 
tiable. Moreover, its derivative is 

dTt , Tt+h — Tt , Tfj — 1 ^ ^ ^ 
— = hm = lim — Tt = LTf 

dt h^o h h^o h 

Q.E.D. 
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Theorem 2.3. The exponential operator T{t) = e^* is the only solution of 
the differential problem 

^ - LT{t), t e R+, 



dt 

T{0) = 1. 

Proof. From the definition of exponential given by Eq. (|3]), it is clear that 
T(0) = 1, and because of the commutativity of the exponents we find 



lim — = ( lim - — 1 e^*. 



dt h^o h \h^o h 

By expanding the exponential one obtains 

, e^^ - 1 , t + Lh + O(h^) - 1 

lim ; = lim ; = L. 

h^O h h-^O h 

At this point, it just remains to prove uniqueness. Let us consider another 
function S{t) and assume that it also satisfies the differential problem. Then 
we define 

Q{s) =T{s)S{t- s), fort>s>0, 
for some fixed t > 0, so Q{s) is differentiable with derivative 
dQ{s) 



ds 



LT{s)S{t-s)-T{t)LS{t-s) = LT{s)S{t-s)-LT{t)S{t-s) = 0. 



Q{s) is therefore a constant function. In particular Q{s) = Q{0) for any 
t > s > 0, which implies that 

T(t) = T{t)S{t -t) = Q{t) = g(o) = S{t). 

Q.E.D. 

Corollary 2.1. Any uniformly continuous semigroup can be written in the 
form Tt = T(t) = e^*, where L is called the generator of the semigroup and 
it is the only solution to the differential problem 

To = 1. 

Given the omnipresence of differential equations describing time evolution 
in physics, it is now evident why semigroups are important. If we can extend 
the domain of t to negative values, then Tt forms a one-parameter group 
whose inverses are given by T_t, so that TtT_t = 1 for all t G R. 

There is a special class of semigroups which is important for our proposes, 
and that will be the subject of the next section. 
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2.3.1 Contraction semigroups 

Definition 2.7. A one-parameter semigroup satisfying ||T(|| < 1 for every 
t > 0, is called a contraction semigroup. 

Now one question arises, namely, which properties should L fulfill to 
generate a contraction semigroup? This is an intricate question in general 
that is treated with more detail in functional analysis textbooks such as [H]- 
1^ . For finite dimensional Banach spaces, the proofs of the main theorems 
can be drastically simplified, and they are presented in the following. The 
main condition is based in properties of the resolvent set and the resolvent 
operator. 

Theorem 2.4 (Hille-Yoshida) . A necessary and sufficient condition for a 
linear operator L to generate a contraction semigroup is that 

1. {A, Re A > 0} C q{L). 

2. WRdm < (ReA)-i. 

Proof. Let L be the generator of a contraction semigroup and A a complex 
number. For a > b we define the operator 

b 



J a 



Since the exponential series converges uniformly (to see this one uses again 
tools of the elementary analysis, specifically the M-test of Weierstrass |13], on 
the Banach space ^B) it is possible to integrate term by term in the previous 
definition 

£(a,6)= / Y.'- -^dt = E ^, (/ 'dt 



E 



;-Al + L)" 



n=0 



n + 1 n + 1 



Thus by multiplying C{a, h) by (Al — L) we get 

(Al - L)£(a, h) = C{a, b){Xl - L) = [e^-^^+'^)- - e(-^i+^)^] . 
Now assume Re A > and take the limit a — 0, 6 — > oo, then as 
hm ||e(-"^+^)''|| = hm ||e-""e^i < lim ||e-""|| ||e^''|| 

< lim ||e-^"|| = 0, 
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and 

we obtain 

(Al - L) 



Ala „La 



lim e """e 

a-)-0 



{XI -L) = 1. 



So if ReA > then A G g{L) and from the definition of the resolvent operator 
we conclude 



Moreover, 



\RLiX)\\ 



< 



Ie""^*||dt = 



< 



-Xtt^Lt 



-At 



dt 



-Rc(A)t^^ 



e""^*||||e"||dt 
1 



Re(A)' 



Conversely, suppose that L satisfies the above conditions U\ and [B Let A 
be some real number in the resolvent set of L, which means it is also a real 
positive number. We define the operator 

Lx = XLRl{X), 

and it turns out that La — )■ L as A — ?■ oo. Indeed, it is enough to show that 
Ai?i,(A) — )■ 1. Since (0, oo) C g{L) there is no problem with the resolvent 
operator taking the limit, so by writing XRl{X) = LRl{X) + 1 we obtain 

lim \\XRl{X) - 1|| = lim \\LRl{X)\\ < lim ||L|| ||i?i(A) || < lim ^ = 0, 

A— >oo A— >oo A— ^oo A— >oo A 

where we have used [E Moreover, 



JA2i?i(A)-Al]t| 



and therefore 



is a contraction semigroup. 



lim e^^ 

A— >oo 



Q.E.D. 
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Apart from imposing conditions such that L generates a contraction semi- 
group, this theorem is of vital importance in general operator theory on ar- 
bitrary dimensional Banach spaces. If the conditions of the theorem 12.41 are 
fulfilled then L generates a well-defined semigroup. Note that for general 
operators L the exponential cannot be rigourously defined as a power se- 
ries, and more complicated tools are needed. On the other hand, to apply 
this theorem directly can be complicated in many cases, so one needs more 
manageable equivalent conditions. 

Definition 2.8. A semi-inner product is an operation in which for each pair 
Xi, X2 of elements of a Banach space 23 there is an associated complex number 
[xi, 0:2], such that 

[xi, Xx2 + /ixa] = \[xi,X2] + X3], (7) 
[xi,xi] = (8) 
|[xi,X2]| < ||xi||||x2||. (9) 

for any Xi, X2, X3 G 53 and A, /i G C. 

Definition 2.9. A linear operator L on 53 is called dissipative if 

Re[x,Lx]<0 VxG*B. (10) 

Theorem 2.5 (Lumer-Phillips) . A necessary and sufficient condition for a 
linear operator L on 03 (of finite dimension ) to be the generator of a con- 
traction semigroup is that L be dissipative. 

Proof. Suppose that {Tt,t > 0} is a contraction semigroup, then just by 
using the properties of a semi-inner product we have that: 

Re[x, {Tt — l)x] = Re[x,Tix] — < |[x,Ttx]| — ||a;||^ 

< ||x||||Tfa;|| - ||xf < ||xf ||Ti|| - \\xf < 0, 

as ||T(|| < 1. So L is dissipative 

Re[a;, Lx] = lim j Re[x, {Tt - l)x] < 0. 

Conversely, let L be dissipative. Considering a point within its spectrum 
A G cr(L) and a corresponding eigenvector x G ker(Al — L). Since [x, (Al — 
L)x] = we have 

Re[x,{Xl - L)x] =Re{X)\\xf -Re[x,Lx] = 0, 
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given that Re[a;, Lx\ < 0, Re(A) < 0. Thus any complex number Re(A) > 
is in the resolvent of L. Moreover, let A G q{L), again from the properties of 
a semi-inner product we have that 

Re(A)||x||^ = Re(A)[x,x] < Re(A[x, x] - [x, Lx]) = Re[x, (Al - 

< |[x, (Al-L)x]| < ||x||||(Al -L)x||, 

so II (Al — L)x\\ > Re(A)||x||, Vx 7^ 0. Then, by setting x = (Al — L)x in the 
expression for the norm of Rl{X): 

II D f\\u Pl(A)x|| ||x|| 

||-rLL(Aj|| = sup — = SUp 



xe<Z,x^O xe!B,x^ker{Al-L) ||(A1 — i^)x|| 

^ II5II 1 



Re(A)||x|| Re(A)^ 



and therefore L satisfies the conditions of the theorem 12.41 

Q.E.D. 



2.4 Evolution families 

As we have seen semigroups arise naturally in the context of problems in- 
volving linear differential equations, as exemplified by Eq. (|6]), with time- 
independent generators L. However some situations in physics require solving 
time-inhomogeneous differential problems, 

r f = L(t)x, xG® 

\ x{to) = Xo, 

where L{t) is a time-dependent linear operator. Given that the differential 
equation is linear, its solution must depend linearly on the initial conditions, 
so we can write 

x{t) = T(j,t„)x(to), (11) 

where T(^t,to) is a linear operator often referred to as the evolution operator. 
Obviously, T(^to,to) = 1 aiid the composition of two consecutive evolution 
operators is well defined; since for t2>ti> to 

x(ti) = T(tj,t(,)x(to) and x{t2) = T(^t2,h)x{ti). 
The evolution operator between t2 and to is given by 

^(t2,<o) = T{t2,h)T{ti,to), (12) 
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and the two-parameter family of operators ^(4^5) fulfills the properties 



T(t,s) = T(t,r)T(r,s), t>r>s , . 

These operators are sometimes called evolution families, propagators, or fun- 
damental solutions. One important fact to realize is that if the generator 
is time- independent, L, the evolution family is reduced to a one-parameter 
semigroup in the time difference T(^t,s) = Tr=t-s- However, in contrast to gen- 
eral semigroups, it is not enough (even in finite dimension) to require that 
the map [t, s) ^ T(^t,s) be continuous to ensure that it is differentiable [17], al- 
though in practice we usually assume that the temporal evolution is smooth 
enough and families fulfilling Eq. f lT^ will also considered to be differen- 
tiable. Under these conditions we can formulate the following result. 

Theorem 2.6. A differentiable family T(t,s) obeying the condition / f73]) is the 
only solution to the differential problems 



dt 



L{t)T(^t,s), t > s 



and 

dT^ _ 



^ - -T^t,s)L{t), t > s 

Proof. If the evolution family T(^t,s) is differentiable 

^ = hm = hm ^^^^^f "V ,,) = L(t)T,,). 

dt h^o h h^o h \ / ^ ' ^ 

Analogously one proves that 



ds 



-Tus)L{t). 



To show that this solution is unique, as in theorem 12. 3[ let us assume that 
there exist some other S{t^s) solving both differential problems. Then, for 
every fixed t and s, we define now Q{r) = T(^t,r)S(r.s) for t > r > s. It is 
immediate to show that ^^^M = and therefore 

dr 

T(t,s) = Q{s) = Q{t) = 

Q.E.D. 
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Again opposed to the case of semigroups, to give a general expression for 
an evolution family T(^t,s) in terms of the generator L{t) is complicated. To 
illustrate why let us define first the concept of time-ordering operator. 

Definition 2.10. The time-ordered operator T of a product of two operators 
L{ti)L{t2) is defined by 



TL{t,)L{t2) = 9{ti - t2)L{t,)L{t2) + 6{t2 - t,)L{t2)L{t 



1^ 



where 9{x) is the Heaviside step function. Similarly, for a product of three 
operators L{ti)L{t2)L{t^), we have 



rL{t,)L{t2)L{ts) = e{t^ - 


- t2)e{t2 


-h] 


)L{t,)L{t2)L{ts) 


+ e{t^ - 


- h)0{h 




)L{t^)L{ts)L{t2) 


+ 9it2 - 


-ti)e{t. 


-ts] 


]L{t2)L{t{)L{t^) 


+ e{t2 - 


- h)e{t^ 


-h] 


]L{t2)L{t^)L{t{) 


+ Q[h - 


-ti)9{ti 


-t2] 


)L{ts)L{t,)L{t2) 


+ Q{t; - 


- t2)e{t2 


-h] 


)L{t^)L{t2)L{t,) 


so that for a product of k operators 


L{tr) ■ ■ 


■L{t, 


;) we can write: 


rL(ti) ■ ■ ■ = ^ - t^(2) 


]---e[t^ 


{k-i) 


— t-K{k)\L\t',T{l)\ ■ 



■ -^[^7r(fc)], 

(14) 

where vr is a permutation of k indexes and the sum extends over all k\ different 
permutations. 

Theorem 2.7 (Dyson Series). Let the generator L{t') to he hounded in the 
interval [t, s], i.e. ||L(t')|| < oo, t' G [t, s]. Then the evolution family Ti^t,s) 
admits the following series representation 

T{t,s) = t + J2 / ■■■/ L{t,)---L{tm)dtm---dti, (15) 

where t > ti > . . . > tm > s, and it can he written symholically as 

T{t, s) = Te/^' = t + J2—i ■■■ '^Lit,) ■ ■ ■ L{tm)dti ■ ■ ■ dt^. 

m=l 

(16) 
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Proof. To prove that expression f[T^ can be formally written as fITBl) note 
that by introducing the time-ordering definition Eq. f lT^ . the k^^ term of 
this series can be written as 



l\ I ■■■ I TL{h)---L{tk)dh---dtk 

= -^2^ / ■■■/ L[t^(^i)] ■ ■ ■ L[t^^k)]dti ■ ■ ■ dtk 



tk-1 

L{ti)---L{tk)dtk---dtu 

where in the last step we have used that t > t7r(i) > • • • > ^7r(fc) > s, and that 
the name of the variables does not matter, so the sum over vr is the addition 
of the same integral k\ times. 

To prove Eq. (IT^ . we need first of all to show that the series is convergent, 
but this is simple because by taking norms in f|T6l) we get the upper bound 



\t — s\ 



ml 

m=0 



exp 



sup \\L{t')\\\t - s\ 

t'e[s,t] 



which is convergent as a result of L(t') being bounded in [t,s]. Finally, the 
fact that (fT5|) is the solution of the differential problem is easy to verify by 
differentiating the series term by term since it converges uniformly and so 
does its derivative (to see this one uses again the M-test of Weierstrass |13] 
on the Banach space ^). Q.E.D. 

This result is the well-known Dyson expansion, which is widely used for 
instance in scattering theory. We have proved its convergence for finite di- 
mensional systems and bounded generator. This is normally difficult to prove 
in the infinite dimensional case, but the expression is still used in this context 
in a formal sense. 

Apart from the Dyson expansion, there is an approximation formula for 
evolution families which will be useful in this work. 

Theorem 2.8 (Time-splitting formula). The evolution family T(^t,s) can he 
given by the formula 



T^t,s) = hm n e(*^+^-*^)^(*^) for t = t^ > . . . > to = s, (17) 

max t,- 1— >0 

j=n-l 

where L[tj) is the generator evaluated in the instantaneous time tj (note the 
descending order in the product symbol). 
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Proof. The idea for this proof is taken from |3S], more rigorous and general 
proofs can be found in [HI H^. Given that the product in f[T7|l is already 
temporally ordered, nothing happens if we introduce the temporal-ordering 
operator 



T(i,,) = lim r TT e(*^+^-*^)^(*^-\ 

' max |— >0 

j=n-l 

— lim 7~eEj'=J(*j+i-*j)-f'(*j) 

max fj |— >0 

where one realizes that the time-ordering operator is then making the work 
of ordering the terms in the exponential series in the bottom term in order to 
obtain the same series expansion as in the top term. Now it only remains to 
recognize the Riemann sum [13] in the exponent to arrive at the expression 

which is Eq. i^^. Q.E.D. 

In analogy to the case of contraction semigroups, one can define contrac- 
tion evolution families. However, since we have now two parameters, there 
is not a unique possible definition, as illustrated below. 

Definition 2.11. An evolution family T(^t,s) is called contractive if ||T(f^s) || < 1 
for every t and every s such that t > s. 

Definition 2.12. An evolution family T(^t,s) is called eventually contractive 
if there exist a Sq such that < 1 for every t > sq. 

It is clear that contractive families are also eventually contractive, but 
in this second case there exist a privileged initial time so, such that one can 
have ||T(t ,,)|| > 1 for some s ^ sq. This difference is relevant and will allow 
us later on to discriminate universal dynamical maps given by differential 
equations which are Markovian from those which are not. 

3 Time evolution in closed quantum systems 

Probably time evolution of physical systems is the main factor in order to 
understand their nature and properties. In classical systems time evolution 
is usually formulated in terms of differential equations (i.e. Euler-Lagrange's 
equations, Hamilton's equations, Liouville equation, etc.), which can present 
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very different characteristics depending on which physical system they cor- 
respond to. From the beginning of the quantum theory, physicists have been 
often trying to translate the methods which were useful in the classical case to 
the quantum one, so was that Erwin Schrodinger obtained the first quantum 
evolution equation in 1926 [i9] . 

This equation, called Schrodinger's equation since then, describes the 
behavior of an isolated or closed quantum system, that is, by definition, a 
system which does not interchanges information (i.e. energy and/or matter 
[50] ) with another system. So if our isolated system is in some pure state 
\ipit)) G "H at time t, where 7i denotes the Hilbert space of the system, 
the time evolution of this state (between two consecutive measurements) is 
according to 

j^\m) = -lH{t)\m), (18) 

where H{t) is the Hamiltonian operator of the system. In the case of mixed 
states p{t) the last equation naturally induces 

^^4|//W.pW1. (19) 

sometimes called von Neumann or Liouville-von Neumann equation, mainly 
in the context of statistical physics. From now on we consider units where 
h=l. 

There are several ways to justify the Schrodinger equation, some of them 
related with heuristic approaches or variational principles, however in all of 
them it is necessary to make some extra assumptions apart from the postu- 
lates about states and observables; for that reason the Schrodinger equation 
can be taken as another postulate. Of course the ultimate reason for its 
validity is that the experimental tests are in accordance with it [^THSB] . In 
that sense, the Schrodinger equation is for quantum mechanics as fundamen- 
tal as the Maxwell's equations for electromagnetism or the Newton's Laws 
for classical mechanics. 

An important property of Schrodinger equation is that it does not change 
the norm of the states, indeed 

= tm)\H\t)\i;it)) - t{m\Hit)m)) = o, (20) 

as the Hamiltonian is self-adjoint, H(t) = H\t). In view of the fact that the 
Schrodinger equation is linear, its solution is given by an evolution family 
f/(t, r), 

m)) = u{tM)mo)), (21) 



19 



where U{tQ, to) = 1- 

The equation fl2Ul) imphes that the evolution operator is an isometry (i.e. 
norm-preserving map), which means it is a unitary operator in case of finite 
dimensional systems. In infinite dimensional systems one has to prove that 
U{t, to) maps the whole "H onto in order to exclude partial isometries, this 
is easy to verify from the composition law of evolution families [51] and then 
for any general Hilbert space 

U\t, to)U{t, to) = U{t, to)U\t, to) = t^ = Ul (22) 

On the other hand, according to these definitions for pure states, the 
evolution of some density matrix p{t) is 

p{ti) = U{tuto)p{to)U\ti,to). (23) 
We can also rewrite this relation as 

piti)=l^ihfyMto)l (24) 
where l^{ti,to) is a linear (unitary ^3]) operator acting as 

U^tuto)[-] = U{tuto)[-]UHtuto). (25) 

Of course the specific form of the evolution operator in terms of the 
Hamiltonian depends on the properties of the Hamiltonian itself. In the sim- 
plest case in which H is time- independent (conservative system), the formal 
solution of the Schodinger equation is straightforwardly obtained as 

[/(t,to) =e-*(*-*°)^/^ (26) 

When H is time-dependent (non-conservative system) we can formally write 
the evolution Dyson expansion, 

f/(t,to) = re^t>(*')'^*'. (27) 

Note now that if H{t) is self-adjoint, —H{t) is also self-adjoint, so in 
this case we have physical inverses for every evolution. In fact when the 
Hamiltonian is time-independent, U{ti,to) = U{ti — to). So the evolu- 
tion operator is not only a one-parameter semigroup, but a one-parameter 
group, with an inverse for every element U^^{t) = U{—t) = (/"^{t), such 
that U~^{t)U{t) = U{t)U~^{t) = 1. Of course in case of time- dependent 
Hamiltonians, every member of the evolution family f/(t, s) has the inverse 
U~^{t,s) = U{s,t) = W{t,s), being this one also unitary. 
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If we look into the equation f l2^ . we can identify p as a member of the 
Banach space of the set of trace-class self-adjoint operators, and U{ti,to) 
some linear operator acting on this Banach space. Then everything con- 
sidered for pure states concerning the form of the evolution operator is also 
applicable for general mixed states, actually the general solution to the equa- 
tion din]) can be written as 

p{t) = W(t,to)p(to) = re^4 ^'"^^ p(to), (28) 

where the generator is given by the commutator Ct{-) = —i[H{t), ■]. This is 
sometimes called the Liouvillian, and by taking the derivative one immedi- 
ately arrives to (|T9l) . 

4 Time evolution in open quantum systems 

Consider now the case where the total Hilbert space is decomposed into two 
parts, in the form "H = Ha '^'Hb, each subspace corresponding to a certain 
quantum system; the question to answer is then, what properties does the 
time evolution of each subsystem have?. 

First of all, one must establish the relation between the total density 
matrix p & H and the density matrix of a subsystem, say A, which is given 
by the partial trace over the other subsystem B: 

PA = Trs(p). (29) 

There are several ways to justify this; see for example [381E6]- Let us focus 
first in the description of a measurement of the system T-La viewed from the 
whole system = Ha ^H-b- For the system T-La a measurement is given by 
a set of self-adjoint positive operators {M^} each of them associated with a 
possible result k. Suppose that we are blind to the system B, and we see the 
state of system A as some pA- Then the probability to get the result k when 
we make the measurement {M^} will be 

PA{k) = TiiMkPA). 

Now one argues that if we are viewing only the part A through the measure- 
ment Mfc, we are in fact observers of an extended system "H through some 
measurement such that, for physical consistency, for any state of the 
composed system p compatible with pA which is seen from A (see figure [1]), 
the following holds 

PA{k) = Tr (Mfcp). 
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But assuming this is true for every state p^, we assume that is a genuine 
measurement operator on "H, which imphcitly imphes that is independent of 
the state p of the whole system on which we are making the measurement. 
Now consider the special case where the state of the whole system has the 
form p = pA®pBi then the above equation is clearly satisfied with the choice 

PA{k) = Tr(Mfcp^) = Tr[(Mfc ® l)p^ ® pe] = Tr(MfcpA) Tt{pb) = Tr(MfcpA). 

Actually this is the only possible choice such that Mk is independent of the 
particular states pa and pb (and so of p). Indeed, imagine that there exist 
another solution, M^, independent of p, which is also a genuine measure on 
'H, then from the linearity of the trace we would get 

PA{k) = TiiMkPA) = Tr(Mfc ® Ip) = Tr(Affcp) Tr[(Mfc - ® l)p] = 0, 

for all Pa and so for all p. This means 

Tr[(Mfc -Mk® l)cj] = 0, 

for any self-adjoint operator a, therefore as the set of self-adjoint operators 
forms a Banach space this implies Mk — Mk ® 1 = and then Mk = Mk ® 1. 



<I> 



A 




B 



<^ <l> 



A 



Figure 1: On the left side we perform a measure over the combined system A 
and B, which we denote by Mk, but we have "disconnected" the information 
provided hj B. If one assume that this situation is equivalent to a measure 
Mk only the system A, this is illustrated on the right side, then Mk = Mk ® 1 
(see discussion on the text). 

Once established this relation, one can wonder about the connection be- 
tween Pa and p. If {\ipf)\ipi)} is a basis of "H then 

PA{k) = Tr[{Mk®t)p] = J2{ijf\{^fmk®l)p\^fM) 

= J2{i;f\Mk{i;f\p\^fMf) =TriMkPA), 
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hence pA is univocally given by the partial trace 

PA = 5^(V^f|p|^f)=TrB(p). 

One important property of the partial trace operation is that, even if 
p = \4'){ip\ was pure, pA can be mixed, this occurs if \ip) is entangled. 

Given the total density matrix pito), by taking the partial trace in (123|) . 
the state of A at time ti, Pa(^i); is given by 

PAih) = TTB[U{ti,to)p{to)U^{h,to)]. (30) 

Provided that the total evolution is not factorized ?7(ti, to) = UA(ti, to) ® 
UB(ti,tQ), both quantum subsystems A and B are interchanging information 
with each other (they are interacting), and so they are open quantum systems. 



4.1 Dynamical maps 

Let us set A as our system to study. A natural task arising from the equation 
(130|) is trying to rewrite it as a dynamical map acting on Ha which connects 
the states of the subsystem A at times to and ti. 

%,to) : pAito) pAih). 

The problem with this map is that in general it does not only depend on the 
global evolution operator U (ti, to) ^'^^ on the properties of the subsystem B, 
but also on the system A itself. 

In order to clarify this, let us write the total initial state as the sum of 
two contributions [571: 



p(to) = PAito) (g) PBito) + Pcorrito), (31) 

where the term pcorrito) is not a quantum state; satisfies 

TrA[Pcor-r(to)] = Tr^ [p^orr (to)] = 0, 

and contains all correlations (classical as well as quantum) between the two 
subsystems. The substitution of (l3Ti) in (l30l) provides 

PA{ti) = TTB{U{tuto)[pA{to)®pB{to)+Pcorr{to)]U\tuto)} 

= TTB{U{ti,to)[pA{to) ® m]U^t,,to)} 

i 

+ TTB[Uiti,to)Pcorrito)U^ih,to)] 

= J2 Ka{tuto)pA{to)Kl{tuto) + 5p(tl, to) 

a 

^ %,,„)[pA(to)], (32) 
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where a = {i,j}, A'{jj}(ti, to) = y/K{'ipj\U(ti,to)\ipi) and we have used 
the spectral decomposition of pBito) = Xli Apparently the op- 
erators Ka(ti,to) depend only on the global unitary evolution and on the 
initial state of the subsystem B, but the inhomogeneous part 6p(ti,to) = 
TTB[U{ti,to)pcorr(to)U^ (tiyto)] niay no longer be independent of pA because 
of the correlation term pcorr [SH]- That is, because of the positivity require- 
ment of the total density matrix p(to); given some pyi(to) and Pb(^o) ^ot any 
Pcorrito) is allowed such that p(to) given by fl3T]) is a positive operator; in the 
same way, given pa(^o) and Pcorr{to) not any psito) is allowed. 

A paradigmatic example of this is the case of reduced pure states. Then 
the following theorem is fulfilled. 

Theorem 4.1. Let one of the reduced states of a bipartite quantum system 
be pure, say pA = \'>P)a{'>P\, then pcorr = 0. 

Proof. Indeed, let us consider first that p = \ip){ip\ is pure as well. By 
using the Schmidt decomposition of \ip) = Aj|uj, i^j), where for all i in 
the sum Aj > 0, and and {|fj)} are orthonormal basis of Ha and Hb 

respectively, we obtain 



Pa = Tri 



Y,^^K){u.\. (33) 



It is elementary to show that a pure state cannot be written as non-trivial 
convex combination of another two states [591160] : then one coefficient, say 
Aj' must be equal to one, and Aj = for all i ^ i'. This implies that p has 
the product form 



p = \Ui>){Ui> \ ® \Vi'){Vi, 

If p is mixed, it can be written as convex combination of pure states 

P = X]Pfc|^fc)(^fc|, 

k 



if p]^ = TTA{\ipk){'^k\), the state of A is 



Pa 



but as Pa is pure, we have only two possibilities, either some p^r = 1 and 
Pk = for k ^ k' or p^^ is the same state pA for every k. In the first 
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case p is pure obviously, for the second option let us introduce the Schmidt 
decomposition again for each state \il)k) = Xli "^j.fc)- 

Pa = P^A = TrA(|V^fc)(^fe|) = ^Xlk\ui,k){ui,k\^ V/c 

i 

SO both k and \ui±) can be taken to be independent of k, and 

Pa = ^ A^|Mi)(Mi|. 

i 

This equation is actually the same as (1331) and so one Aj' is equal to one, and 
Aj = for all i 7^ i'. With these requirements we get \ipk) = l^i', fj'.fc), and 
finally 

P = ^Pk\i^k){i^k\ = y^^Pk\ui',Vi',k){ui',Vi',k\ 

k k 

= \ui'){ui' \ ® ^Pk\vi',k){vi',k\ = \tp)A{i^\ ® Pb, 

k 

and hence pcorr = 0. 

Q.E.D. 

As a result, the positivity requirement forces each term in the dynamical 
map to be interconnected and dependent on the state upon which they act; 
this means that a dynamical map with some values of /^^(ti, to) and 5p(ti, to) 
may describe a physical evolution for some states pyi(to) and an unphysical 
evolution for others. 

Given this fact, a possible approach to the dynamics of open quantum 
systems is by means of the study of dynamical maps and their positivity 
domains [GTHM], however that will not be our aim here. 

On the other hand, the following interesting result given independently 
by D. Salgado et al. [HZ] and D. M. Tong et al. [BB] provides an equivalent 
parametrization of dynamical maps. 

Theorem 4.2. Any kind of time evolution of a quantum state pA can always 
be written in the form 

PAih) = '^Ka (ti, to, Pa) pA{to)Kl (ti, to, pa) , (34) 

a 

where iti,to,pA) are operators which depend on the state pa at time to- 



25 



Proof. We give here a very simple proof based on an idea presented in 
Let Pa(^o) and p^(ti) be the states of a quantum system at times to and ti 
respectively. Consider formally the product ^^(^0) 8) G "Ha ®'Ha, and 
a unitary operation f/swAP which interchanges the states in a tensor product 

t^SWApl^l) ® 1^2) = |V^2) ® 1^1), SO ;7sWAp(^ ® ^)^SWAP = B ® A. It is 

evident that the time evolution between ^^(^0) and pAiti) can be written as 



^{ti,to)[pA{to)\ = Tt2 UswAPpA{to) ® PA(il)f/swAP = P^(* 



where Tt2 denotes the partial trace with respect to the second member of 
the composed state. By taking the spectral decomposition of Pyi(ti) in the 
central term of the above equation we obtain an expression of the form (l34l) . 

Q.E.D. 

Note that the decomposition of 0341) is clearly not unique. For more details 
see references [67] and [68] . 

It will be convenient for us to understand the dependence on pA of the 
operators Ka (ti, to, pa) in flMl) as a result of limiting the action of fl32l) . That 
is, the above theorem asserts that restricting the action of any dynamical map 
fl5^ from its positivity domain to a small enough set of states, actually to 
one state "p^i" if necessary, is equivalent to another dynamical map without 
inhomogeneous term Sp{ti,to). However this fact can also be visualized as 
a by-product of the specific nonlinear features of a dynamical map as (132|) 

[SHIED]. 

As one can already guess, the absence of the inhomogeneous part will 
play a key role in the dynamics. Before we move on to examine this in the 
next section, it is worth remarking that C. A. Rodriguez et al. [71] have 
found that for some classes of initial separable states p(to) (the ones which 
have vanishing quantum discord), the dynamical map fl3ll) is the same for a 
set of commuting states of pa- So we can remove the dependence on pA of 
Ka iti,to, Pa) when applying flM|l inside of the set of states commuting with 

PA- 



4.2 Universal dynamical maps 

We shall call a universal dynamical map (UDM) to a dynamical map which is 
independent of the state it acts upon, and so it describes a plausible physical 
evolution for any state pa- In particular, since for pure pA(to), Pcorrito) = 0, 
according to (152]) the most general form of a UDM is given by 

PA(ti) = £(t,,to)[pAito)] = ^Ka{ti,to) pA{to)Kl{ti,to) , (35) 

a 
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in addition, the equality 



Y,Kiiti,to)K^{tuto) = t (36) 

a 

is fulfilled because Tr[p^(ti)] = 1 for any initial state pa(^o)- Compare ( l34|) 
and ([35]). 

An important result is the following (see figure |2]): 

Theorem 4.3. A dynamical map is a UDM if and only if it is induced from 
an extended system with the initial condition p(to) = Pa(^o) ® Psito) where 
psito) is fixed for any pA{to). 

Proof. Given such an initial condition the result is elementary to check. 
Conversely, given a UDM fl35|) one can always consider the extended state 
p(to) = pA{to) ® \^Jj)b{^\, with \4')b G 'Hb fixed for every state pyi(io), and 
so Ka(ti,to) = B{(pa\U{ti,to)\ip) B, whcrc \(f)a)B IS & basis of Hb- Since 
-f^a (^i5^o) only fixes a few partial elements of U{ti,to) and the requirement 
fl36|l is consistent with the unitarity condition 

J2B{^\U\ti,to)\<j)o.)B{<Pa \Uiti,to)\4')B = B{'ip\U'' (ti,to)U(ti,to)\lp) B 
a 

= B{i^\t(^t\i^)B = t, 

such a unitary operator f/(ti,to) always exists. 

Q.E.D. 



U{ti,to) 

PA{to) ® pBiU) — ^ pAB{tl) 



Pa(*o) 



'(ti,*o) 



PA{tl) 



Figure 2: Schematic illustration of the result presented in theorem 14.31 

Moreover UDMs possess other interesting properties. For instance they 
are linear maps, that is for pA{to) = -^jPaH^o) (with Aj > and Aj = 1) 
a straightforward consequence from definition of UDM is 
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In addition, UDMs are positive operations (i.e. the image of a positive op- 
erator is also a positive operator). 

Furthermore, a UDM holds another remarkable mathematical property. 
Consider the following "Gedankenexperiment" ; imagine that the system AB 
is actually part of another larger system ABW, where the W component 
interact neither with A nor with B, and so it is hidden from our eyes, playing 
the role of a simple bystander of the motion of A and B. Then let us consider 
the global dynamics of the three subsystems. Since W is disconnected from A 
and B, the unitary evolution of the whole system is of the form UABiti, to) ® 
Uwitiyto), where Uwiti,to) is the unitary evolution of the subsystem W. If 
the initial state is PABwito) then 

PABwiti) = UAB{ti,to) ® Uw(ti,to)pABw{to)U\^(ti,to) (g) f/^(tl,to)- 

Note that by taking the partial trace with respect to W we get 

PAB(ti) = UABiti, to)TTw Uw{ti,to)pABw{to)Uly(ti,to) Ul^^(ti,to) 

= UAB(ti,to)TTw Uly(ti,to)Uw(ti,to)pABw(to) ^i_B(^l!^o) 
= UAB(ti,to)pAB{to)U\^(ti,to), 

so indeed the presence of W does not perturb the dynamics of A and B. 
Since we assume that the evolution of A is described by a UDM, according 
to the theorem 14.31 the initial condition must be PAsito) = Pa(^o) ® PB{to), 
and then the most general initial condition of the global system which is 
compatible with that has to be of the form pABwito) = Paw (to) (8 PB(to), 
where ^^(to) is fixed for any pAwito). Now let us address our attention to 
the dynamics of the subsystem AW: 

pAwih) = TrB [UABiti,to) ® Uwih, to)pAw{'to) ® pi?(to) 

by using the spectral decomposition of PB{to), as we have already done a 
couple of times, one obtains 

pAwih) = ^Kc,{ti,to) ®Uw{ti,tQ)pAw{tQ)Kl{ti,tQ) (^U'^iti.tQ) 

a 

= £{ti,to) ® ^itiM) [PAwiio]] ■ 

Here £(t-i_,to) is a dynamical map on the subsystem A and the operator lA[ti,to) ["] = 
f/vK(^i5 ^o)[']^vi/(^i! ^o) is the unitary evolution of W. As we see, the dynami- 
cal map on the system AW is the tensor product of each individual dynamics. 
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So if the evolution of a system is given by a UDM, £^{ti,to), then for any uni- 
tary evolution W(ti,fo) of any dimension, S(^ti,to) ®^ii,to) ^^^^ ^ UDM; and 
in particular is a positive-preserving operation. By writing 

as 1 (g> U{t^,to) is a unitary operator, the condition for positivity-preserving 
asserts that £{t^,to) ® 1 is positive. Linear maps fulfilling this property are 
called completely positive maps [72] . 

Since we can never dismiss the possible existence of some extra system 
W , out of our control, a UDM must always be completely positive. This may 
be seen directly from the decomposition f l5^ . The reverse statement, i.e. any 
completely positive map can be written as ( 135|) . is a famous result proved 
by Kraus [7]; thus the formula (135|) is referred as Kraus decomposition of a 
completely positive map. Therefore another definition of a UDM is a (trace- 
preserving) linear map which is completely positive, see for instance [2T|[38]. 

Complete positivity is a stronger requirement than positivity, the best 
known example of a map which is positive but not completely positive is the 
transposition map (see for example [38]), so it cannot represent a UDM. It 
is appropriate to stress that, despite the crystal-clear reasoning about the 
necessity of complete positivity in the evolution of a quantum system, it is 
only required if the evolution is given by a UDM, a general dynamical map 
may not possess this property [CTHBH I ^ I ^ ITU]. Indeed, if 5p(ti,to) 7^ 0, 
^{ti,to) given by ( 132|) is not positive for every state, so certainly it cannot be 
completely positive. In addition note that if pABwiio) 7^ PAwiio) ® PB(to) 
we can no longer express the dynamical map on AW as i^(ti,fo) ®^ti,to)' 
the physical requirement for complete positivity is missed. 

4.3 Universal dynamical maps as contractions 

Let us consider now the Banach space 03 of the trace-class operators with 
the trace norm (see comments to the definition 12. ip . Then, as we know, an 
operator p G ^ is a quantum state if it is positive-semidefinite, this is 

and it has trace one, which is equal to the trace norm because the positivity 
of its spectrum, 

p G ^ is a quantum state p > and ||p||i = L 

We denote the set of quantum states by ^+ C ^, actually the Banach space 
03 the smallest linear space which contains ^B"*"; in other words, any other 
Unear space containing also contains 03. 
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We now focus our attention in possible linear transformations on the 
Banach space S{^) : t— ?■ 03, and in particular, as we have said in section 
I2.H the set of these transformations is the dual Banach space 03*, in this 
case with the induced norm 



|£^||i = sup 



l^(p)lh 



pe<s,p^o L llPlli 

Particularly we are interested in linear maps which leave invariant the 
subset i.e. they connect any quantum state with another quantum 
state, as a UDM does. The following theorem holds 

Theorem 4.4. A linear map S G ^* leaves invariant if and only if it 
preserves the trace and is a contraction on 03 

ll^lli < 1- 

Proof. Assume that S G leaves invariant the set of quantum states S"*", 
then it is obvious that it preserves the trace; moreover, it implies that for 
every positive p the trace norm is preserved, so = Let cr € 03 

but cr ^ ^B"*", then by the spectral decomposition we can split a = a+ — a", 
where 

V+ = EA,#,)(^,|, for A, >0, 
= -E Aj#i)(V^j|, for Aj < 0, 

here Xj are the eigenvalues and the corresponding eigenvectors. Note 

that cr"'" and a~ are both positive operators. Because \ are orthogonal 

projections we get (in other words as ||cr||i = J2j 

\W\\i = + 

but then 

||^(a)||i = < + = + = \\a\U, 

so £" is a contraction. 

Conversely, assume that £^ is a contraction and preserves the trace, then 
for p G ^"'" we have the next chain of inequalities: 

||p||i = Tr(p) = Tr[£(p)]<||^(p)||i<||p||i, 

so \\S{p)\\i = IIpIIi and Tr[f (p)] = ||^(p)||i. Since p G ^+ if and only if 
IIpIIi = Tr(p) = 1, the last equality implies that £{p) G for any p G ^"'". 

Q.E.D. 

This theorem was first proven by A. Kosakowski in the references [731171] , 
later on, M. B. Ruskai also analyzed the necessary condition in |75j . 
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4.4 The inverse of a universal dynamical map 

Let us come back to the scenario of dynamical maps. Given a UDM, S(^ti,to)j 
one can wonder whether there is another UDM describing the evolution in 
the time reversed sense £{to,ti) in such a way that 

£{to,ti)£iti,to) = £{ttto)^{ti,to) = 1- 

Note, first of all, that it is quite easy to write a UDM which is not bijective 
(for instance, take {!'?/'-«)} to be an orthonormal basis of the space, then define: 
^(p) = Zln \(P){i'n\p\i^n){4>\ = Tr(p)|0)(0|, being |0) an arbitrary fixed state), 
and so the inverse does not always exist. But even if it is mathematically 
possible to invert a map, we are going to see that in general the inverse is 
not another UDM. 

Theorem 4.5. A UDM, £{ti,to), can be inverted by another UDM if and only 
if it IS unitary £{ti,to) = ^ti.to)- 

Proof. The proof is easy (assuming Wigner's theorem) but lengthy. Let us 
remove the temporal references to make the notation less cumbersome, so we 
write S = S(^ti,to)- Since we have shown that UDMs are contractions on ®, 

\\£icr)\\i<\\a\\i, a €03. 

Let us suppose that there exists a UDM, S~^, which is the inverse of S; then 
by applying the above inequality twice 

||a||i = ||^-i^(a)||i<||^(a)||i<||a||i, 

and we conclude that 

\\Sia)h = \\a\\i, Vae«. (37) 

Next, note that such an invertible S would connect pure stares with pure 
states. Indeed, suppose that \ip){ip\ is a pure state, and £{\ip){ip\) is not, 
then we can write 

£{\ij){i!\)=ppi + {l~p)p2, l>p>0, 

for some pi 7^ P2 7^ But after taking the inverse this equation 

yields 

and is a (bijective) UDM, S~^{pi) 7^ £~^{p2) are quantum states; so such 
a decomposition is impossible as is pure by hypothesis; and a pure 
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state cannot be decomposed as a non-trivial convex combination of another 
two states. 

Now consider a = — \ip2){4'2\) for two arbitrary pure states 

\xl:i){ipi\ and |'?/^2)('?A2|- We can calculate its eigenvalues and eigenvectors 
by restring the computation to the two dimensional subspace spanned by 
and \ip2)', because any other state can be decomposed as a linear 
combination of\ipi), \ip2) and \4'±), such that cr\ip±) =0. Thus the eigenvalue 
equation reads 

a(a|V^i) + /3K'2)) = A(a|^i)+/3|V^2)), 

for some A, a and /3. By taking the scalar product with respect to and 
{1^21, and solving the two simultaneous equations one finds 

so the trace norm ||a"||i = J2j l^jl turns out to be 

||a||i = v/l-|(V'2|V^i)P. (38) 
However, since S connects pure states with pure states, 

and therefore from equations fl37|) and fl38l) we conclude that 

|(V^2|^l)| = |(^2|^l)|. 

Invoking Wigner's theorem [5ll[59], the transformation S have to be imple- 
mented by 

£ = VpV\ 

where is a unitary or anti-unitary operator. Finally since the effect of an 
anti-unitary operator involves the complex conjugation [M], that is equivalent 
to the transposition of the elements of a trace-class self-adjoint operator, 
which is not a completely positive map; we conclude that V has to be unitary, 

Q.E.D. 

From another perspective, the connection between the failure of an inverse for 
a UDM to exist just denotes the irreversibility of the universal open quantum 
systems dynamics. 
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4.5 Temporal continuity. Markovian evolutions 

One important problem of the dynamical maps given by UDM is related with 
their continuity in time. This can be formulated as follows: assume that we 
know the evolution of a system between to and ti, S(^ti,to)j between ti and 
a posterior time t2, £{t2M)- From the continuity of time, one could expect that 
the evolution between Iq and t2 was given by the composition of applications, 
^{t2,to) ~ ^{t2,ti)^{ti,to)- However, the exact meaning of this equation has to 
be analyzed carefully. 

Indeed, let the total initial state be p(to) = PA{to) ® PB(to), and the uni- 
tary evolution operators U{t2,ti), U(ti,to) and U{t2,to) = U(t2,ti)U(ti,to). 
Then the evolution of the subsystem A between to and any subsequent t is 
clearly given by, 

= J2^c.{tuto) PA{to)Kl{tuto) , 

a 

and 

^{h,to)[pA{to)] = Tr B[U{t2,to)pAito) ^ pB{to)U^{t2,to)] 

= Y.K^{t2,t^) PAitQ)Kl{t2M) ■ 
a 

The problem arises with £(t2,ti)'- 

£{t2M)[pA{tl)] = TlB[U{t2M)p{tl)U\t2M)] 

= ^ K„ (ta, Pa) pA{ti)Kl (ta, ti, Pa) , 

a 

being p(ti) = f/(ti,to)[PA(^o) ® PB{to)]U'' (ti,to)- Since p(ti) is not a tensor 
product in general, and it depends on what initial states were taken for both 
subsystems, £{t2M) not the form of an UDM (see figure [3j). 

Let us to illustrate this situation from another point of view. If the 
application £{t'i,to) is bijective it seems that we can overcome the problem by 
defining 

= ^(t2M)^{tlMV ^^^^ 
which trivially satisfies the relation £{t2,to) = £{t2,ti)£{ti,to)- However, as we 
have just seen, is not completely positive unless it is unitary, so in 

general the evolution £{t2,h) is still not a UDM [75] . 

The same idea can also be used to define dynamical maps beyond UDM 
[77] . That is, given the global system in an arbitrary general state p{tQ), 



33 




not UDM 

Figure 3: General situation describing the dynamics of open quantum sys- 
tems in terms of a UDM. At time to the total state is assumed to be factor- 
ized, p{to) = pAiio) ® Psiio), and the dynamical maps from this point are 
UDMs. However it is not possible to define a UDM from some intermediate 
time tl > to; because generally at this time the total state can no longer be 
written as a tensor product. 



if this is not a tensor product, we know that at least in some past time 
tinitiai, before the subsystems started interacting, the state was factorized 
p(iinitiai) = Pyi(^initiai) ® (^initial) • If we manage to find a unitary operation 

t^"^(io,^initial) which maps p(to) to the state p(tinitial) = PA(^initial) ®PB(iinitial) 

for every reduced state pA(to) compatible with p(to)) having pB(tinitiai) the 
same for all pA(to), we can describe any posterior evolution from to to ti as 
the composition of the inverse map and the UDM from tinitiai to ti as in fl39l) . 
However to find such an unitary operator can be complicated, and it will not 
always exist. 

On the other hand, as a difference with the dynamics of closed systems, 
it is clear that these difficulties arising from the continuity of time for UDMs 
make it impossible to formulate the general dynamics of open quantum sys- 
tems by means of differential equations which generate contractive families 
on !B (i.e. families of UDMs). However, sometimes one can write down a 
differential form for the evolution from a fixed time origin to, which is only a 
valid UDM to describe the evolution from that time to. In other words, the 
differential equation generates an eventually contractive family from to on 03 
(see definition 12.121) . 

To illustrate this, consider a bijective UDM, £{t,tQ)i then by differentiation 
(without the aim to be very rigorous here) 

^ = = %^f,7,;„,ip.wi = A[p.(i)i, (40) 

where Ct = This is the ultimate idea behind the "time convolu- 

tionless" method which will be briefiy described in section 17.21 If the initial 
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condition is given at Iq, then, by construction, the solution to the above 
equation is pAit) = S(^t,to)[PA(to)]- Nevertheless, there is not certainty that 
the evolution from any other initial time is a UDM. Indeed, as we may guess, 
it is not difficult to write solutions given the initial condition p{ti) {ti > to) 
as ([39]). 

Definition 4.1. We will say that a quantum systems undergoes a Markovian 
evolution if it is described by a contractive evolution family on ?B and thus 
we recover the composition law for the UDMs 

^(i2,io) = ^{t2,tl)^itl,to), (41) 

which is the quantum analogue to the Chapman-Kolmogorov equation in 
classical Markovian process (see next section [37T]) . This property is sometimes 
also referred as "divisibility condition". 

Typically, the dynamics of an open quantum system is not Markovian, 
because as we have already said both systems A and B develop correlations 
Pcorr during their evolution and so the reduced dynamical map is not a UDM 
for some intermediate starting time ti. However if the term which goes with 
Pcorr does uot affect so much the dynamics, a Markovian model can be a 
good approximate description of time evolution. We will analyze this deeply 
in section O 

5 Quantum Markov process: mathematical 
structure 

Before studying under which conditions can an open quantum system be ap- 
proximately described by a Markovian evolution, we examine in this section 
the structure and properties of quantum Markovian process. 

5.1 Classical Markovian processes 

To motivate the adopted definition of a quantum Markov process, let us recall 
the definition of classical Markov processes. More detailed explanations can 
be found in [1^121(1781179] : here we just sketch the most interesting properties 
for our purposes without getting too concern with mathematical rigor. 

On a probability space, a stochastic process is a family of random variables 
{X{t),t e / C R}. Roughly speaking, a stochastic process is specified by 
a random variable X depending on a parameter t which usually represents 
time. Assume that / is a countable set labeled by n, then a stochastic process 
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is a Markov process if the probability that the random variable X takes a 
value Xn at any arbitrary time is conditioned only by the value that 
it took at time tn-i, and does not depend on the values at earlier times. This 
condition is formulated in terms of the conditional probabilities as follows: 

p{Xn,tn\Xn-l,tn-i;...;Xo,to)=p{Xn,tn\Xn-l), '^U, E I . (42) 

This property is sometimes expressed in the following way: a Markov process 
does not have memory of the history of past values of X. They are so-named 
after the Russian mathematician A. Markov. 

Consider now this process on a continuous interval J, then tn-i is time 
infinitesimally close to So for any t > t' just from the definition of 
conditional probability 

p{x, t] x' , t') = p{x, t\x', t')p{x', t'), 

where p{x, t; x' , t') is the joint probability that random variable X takes value 
X at time t and x' at time t' . By integrating with respect to x' we find the 
relation between unconditional probabilities 



p{x,t) = J dx'p{x,t\x' ,t')p{x' ,t'), 

and let us write it as 
p{x,t) = j dx'K{x,t\x',t')p{x',t') with K{x,t\x',t') = p{x,t\x',t'). (43) 

The Markov process is called homogeneous if K{x,t\x' ,t') is only a function 
of the difference between the two time parameters involved K{x,t\x' ,t') = 
Kt-t'{x\x'). 

Next, take the joint probability for any three consecutive times > ^2 > 
ti and apply again the definition of conditional probability twice 

p{x3,t3;X2,t2]Xi,ti) = p{x3,ts\x2,t2;Xi,ti)p{x2,t2]Xi]ti) 

= p{x3, t3\x2, t2] Xi, ti)p{x2, ^2^1, tl)p(xi, ti). 

The Markov condition implies that p{x3,t3\x2,t2; Xi,ti) = p(x3, t3|x2, ^2), 
then by integrating over X2 and dividing both sides by p{xi,ti) one gets 

p{x3,t3\xi,ti) = J dX2p{x3,t3\x2,t2)p{x2,t2\xi,ti), 

which is called Chapman-Kolmogorov equation. Observe that with the nota- 
tion of (gSD this is 

K{x3,t3\xi,ti) = I dX2K{x3,t3\x2,t2)K{x2,t2\xi,ti). (44) 
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Thinking of K{x, t\x' , t') as propagators of the time evolution between t' and 
t (see Eq. equation flH|) express that they form evolution families with 
that composition law. Moreover, since K{x,t\x' ,t') are conditional probabil- 
ities they connect any probability p{x', t') with another probability p{x, t) so 
in this sense they are universal (preserve the positivity of p{x',t') and the 
normalization) . 

We can now draw a clear parallelism with a quantum setting. In the 
definition of quantum Markov process (definition 14. ip . the role played by the 
probabilities p{x, t) is played by the density operators p{t) and the role of 
conditional probabilities p{x, t\x', t') is played by UDMs £{t,t')^ such that fHTjl 
is the quantum analog to the classic Chapman-Kolmogorov equation (jUl) . 



5.2 Quantum Markov evolution as a differential equa- 
tion 

For positive e, consider the difference 

p{t+e)-p{t) = [£:(i+,,o)-£:(t,o)]p(o) = -i]%o)[p(o)] = [S^^t+e,t)-t\p{t), 

provided that the limit e — )■ is well defined (we assume that time evolution 
is smooth enough). We can obtain a linear differential equation for p{t) 
(called master equation) 

Ml , „„. + ^) - ^ ,.„. I^Ct.^-) - ^p(t) , CMt), (45) 

at e->o e e-i>o e 

where by definition the generator of the evolution is 

A = liml^^^±!^. 

If some quantum evolution is Markovian according to our definition 14.11 
we may wonder what is the form of these differential equations. The answer 
is given by the following theorem, which is the main result of this section. 

Theorem 5.1. A differential equation is a Markovian master equation if and 
only if it can he written in the form 

dp{t) 



dt 



-t[H{t),p{t)] + J2^,{t) V,{t)p{t)V^{t) - l{Vi{t)V,{t),p{t)} 

(46) 

where H{t) and Vk(t) are time- dependent operators, with H{t) self-adjoint, 
and 'jkit) > for every k and time t. 
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Proof. Let us show the necessity first (for this we follow a generalization of 
the argument given in [H] and [21] for time-independent generators). If £{t2,ti) 
is a UDM for any t2 > ti it admits a decomposition of the form 



Actually, let {Fj,j = 1, . . . , A^^} be a complete orthonormal basis with re- 
spect to the Hilbert-Schmidt inner product (Fj,Ffc)Hs = Tr(FjFfc) = 6jk, in 
such a way that we chose F^2 = t/\/N and the rest of operators are then 
traceless. Now the expansion of the Kraus operators in this basis yields 

where the elements 



{t2, h) = J2 ^-(^2, tl))HS (^k, K^{t2, tl)) 



US 



form a positive semidefinite matrix c[t2,ti). Indeed, for any A^^-dimensional 
vector V we have that 

{■V,c{t2,ti)\) = ^V*Cjk{t2,ti)Vk = ^ i^Vk{Fk,Kc{t2,ti))^^\ > 0. 

j,k a \ k / 

Since this holds for any t2 > ti, let us take ti = t and t2 = t + e, then the 
generator reads 



iV 



Cj7v2 (t -f e, t ) cjv2 • (t -f e, t ) . 
— ^ F,p H — -pFl 



+ E 



(47) 



Now, we define the time-dependent coefficients ajk{t) by 

cj\[2N2{t + e,t) — N 



aN2]^2{t) 

ajAr2(t) 

ajk{t) 



lim 

e-S>0 e 

lim^^^^^^^±^, , = l,...,iV^-l, 

e->0 e 



Cjfc(t-Fe,t) 



lim ^^^^^ 
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the following operators, 
G{t) 



and 



= 5^ [^'W - -f (*)] . 



which is self-adjoint. In terms of these operators the generator fj47|) can be 
written like 

= -2[//(t),p] + {G(t),p} + (^At)F3pFl 

j,k=i 

Since the UDMs preserve trace for every density matrix p, 
= Tr[A(p)] = Tr 



2G{t) + J2 <^At)FlF, 
j,k=i 



and therefore we conclude that 



Git) = E ^J>^(t)FlF^ 



j,k=i 



Note that this is nothing but the condition '^^Kl^(t2,ti)Ka(t2,ti) = 1 ex- 
pressed in differential way. Thus the generator adopts the form 



Ct{p) = -t[H{t),p] + J2 ^At) F.pFl - \{FlF,,p] 

j,k=i L 

Finally note that because of the positive semidefiniteness of Cjk{t + e,t), the 
matrix ajk{t), j, k = 1, . . . , N"^ — 1 is also positive semidefinite, so for any t 
it can be diagonalized by mean of a unitary matrix u(t); that is, 

'Yumj(t)ajk(t)ulf,(t) = Jm(t)Smn, 
3,k 

where each eigenvalue is positive 7m (^) > 0. Introducing a new set of opera- 
tors Vfc(t), 

Vkit) = J2 F, = n,,(t)\4(t), 

j=i k=i 
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we obtain the desired result 

A(p) = -z[iJ(t),p] + ^7fcW 



Vk{t)pVi{t)-^{V^it)Vk{t),p} 



with 7fc(t) > 0,\/k,t. 

Assume now that some dynamics is given by a differential equation as 
fj^6|) . Since it is a first order linear equation, there exists a continuous family 
of propagators satisfying 

for every time t2 > ti > t^. It is clear that these propagators are trace pre- 
serving (as for any p, 1i{Ctp) =0). So we have to prove that for any times 
t2 > ti the dynamical maps £{t2,ti) are UDMs, or equivalently completely pos- 
itive maps. For that we use the approximation formulas introduced in section 
[21 firstly the time-splitting formula (fT7|) allow us to write the propagator like 



(t2,ii) 



lim 



n 

j=n-l 



e 



(48) 



where t2 = t'^> t'^_i > ■ ■ ■ > tQ = ti. Next, take any element of this product 



e , Tt = - > 0, 



and since the constants 7A;(t) are positive for every time, we can write this 
generator evaluated at the instant t'^ as 



k 

where Vfc(4) = \/lk{t'i)yk{t'()- Now we split this generator in several parts, 

M 



where 



A;=0 



M 



1 r.~.. 



-^[H{t',),{.)]-Y,7^{vM)vuai-)} 



A(-) = v,{t',){-)v;^{t',), k 



,M, 



(49) 
(50) 
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and here M is the upper hmit in the sum over k in fH^ . On one hand, 
differentiating with respect to r one can easily check that 



which has the form OpO\ so it is completely positive for all r. On the other 
hand, for r > 



m=0 



has also the Kraus form '^^OmpO}^, so it is completely positive for all k. 
Now everything is ready to invoke the Lie- Trotter product formula since 



e I 



= lim TT ( 

TJ. — Von \ J- J- 



is the composition (infinite product) of completely positive maps e'''*'^/", it 
is completely positive for all r. Finally, because this is true for every instan- 
taneous t'l^ the "time-splitting" formula (l48l) asserts that £{t2,h) is completely 
positive for any t2 and ii, because it is just another composition of completely 
positive maps [50] . 

Q.E.D. 

This result can be taken as a global consequence of the works of A. Kos- 
sakowski [73l[711[8T] and G. Lindblad [82], who analyzed this problem for the 
case of time-homogeneous equations, that is, when the UDMs depend only 
on the difference t = t2 — ti, £{t2,ti) = £t and so they form a one-parameter 
semigroup (not a group because the inverses are generally not a UDM, see 
theorem l4.5p . As we know, under the continuity assumption, £r = e^'^ can be 
defined and the generator C is time- independent, moreover it has the form 



Cpit) = -t[H,pit)] + J2lk IvkPmi - l{V^V,,pit)} 



(51) 



This is the result proven by Gorini, Kossakowski and Sudarshan [81] for finite 
dimensional semigroups and by Lindblad [H2| for infinite dimensional systems 
with bounded generators, that is, for uniformly continuous semigroups 



5.3 Kossakowski conditions 

It is probably worthwhile to dedicate a small section to revise the original 
idea of Kossakowski's seminal work [73l[711[8T] , as it has clear analogies with 
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some results in classical Markov processes. The route that A. Kossakowski 
and collaborators followed to arrive at f lHT]) was to derive some analogue to 
the classical Markovian conditions on the generator C That is, for classical 
Markov process in finite dimensional systems, a matrix Q is the generator 
of an evolution stochastic matrix e^"^ for any r > if only if the following 
conditions in its elements are satisfied: 

Qii < 0, Vz, (52) 
Q^, > 0, Vz ^ J, (53) 
^g,, =0, Vz. (54) 



The proof can be found in [75] . To define similar conditions in the quantum 
case, one uses the Lumer-Phillips theorem I2.5[ as Er = e^'^ has to be a con- 
traction semigroup in the Banach space of self-adjoint matrices with respect 
to the trace norm. For that aim, consider the next definition. 

Definition 5.1. Let cr G 03, with spectral decomposition a = Ylj'^jPj- 
define the linear operator sign of a as 

sgn(a) = ^sgn{aj)Pj, 

j 

where sgn((Tj) is the sign function of each eigenvalue, that is 

1, aj > 
sgn(crj) = < 0, = 

-1, CTj < 0. 

Definition 5.2. We define the following product of two elements a, p G 53, 

p]k = \W\\ T¥[sgn(cr)p]. (55) 

Proposition 5.1. The product ( f53]) is a semi-inner product. In addition it 
is a real number. 

Proof. The proof is quite straightforward. It is evident that for the product 
(155|) . ([7]) and ([8]) hold. Finally ([9]) also follows easily. By introducing the 
spectral decomposition of o" = J2j ^j^j i 



W,p]k\ = \W\ 



< ||a||^|sgn(a,)||Tr(i^»| 

j 

< Mj2mp^p)\; 
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and we also use the spectral decomposition of p = J^jPj^jy obtain 



k,j k k 

The reality of [a, p\k is clear from the fact that both sgn(cr) and p are self- 
adjoint. 

Q.E.D. 

Now we can formulate the equivalent to the conditions fl5^ . fl55]l and in 
the quantum case. 

Theorem 5.2 (Kossakowski Conditions). A linear operator C on !B (of finite 
dimension N) is the generator of a trace preserving contraction semigroup 
if and only if for every resolution of the identity = (Pi, P2, ■ ■ ■ , Pn), 1 = 
^jPj, the relations 



Aii{^) < = l,2,...,iV), 

> (^^j = l,2,...,iV), 



N 



Y^A.i^) = (j = l,2,...,iV) 



(56) 
(57) 

(58) 



1=1 



are satisfied, where 



A,m = Tr[P,CiP,)]. 



Proof. The condition 0581) is just the trace preserving requirement Tr[£(o')] 
0, a G 53. We have 



N 



N 



E = E Tr[^.^(^',)] = Tr[£(P,)] = 0, VP, G ^ 

i=l i 

^ Tr[£(a)] = 0, We *B. 
On the other hand, the theorem of Lumer-Phillips 12.51 asserts that 

[a,Cicx)]K<0; (59) 
taking a to be some projector P we obtain 

Tr[P,£(P)] < 0, (60) 
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so flSBj) is necessary. Furthermore (jSH]) can be rewritten as 

k,j j 

+ ^sgnK)a,Tr[P,/:(P,)] <0, (61) 

where we have used the spectral decomposition a = J2j '^jPj again. By 
splitting (jnSD in a similar way 

Tr[P,£(P,)] = -Y, Tr[Pfc/:(P,)], for each j = 1, . . . , iV. 

Here note that on the right hand side the sum is just over the index k (not 
over j). Inserting this result in fl6ip yields 

^afc,Tr[P,X(P,)] >0, (62) 

where the sum is over the two indexes, with 

akj = \(Jj\[l - sgn(crfc) sgn(crj)] > 0. 

Thus for trace-preserving semigroups the condition fl6Tl) is equivalent to fl56|) 
and (1^ . Therefore if fISB]) . flFTI) and flSSl) are satisfied C is the generator of 
a trace preserving contraction semigroup because of fIBUl) and flS2]) . It only 
remains to prove the necessity of (157|) . that can be seen directly from (l62ll . 
However, it is also possible to show it from the definition of generator: 

, , Tr[Pfc(£^ - 1)P-1 , Tr[Pfc£^(P)l 
Tr Pfc£ P, = lim ^ ^ = lim ^ ^ > 0, 

since Sr{Pj) > as £^t- is a trace preserving and completely positive map. 

Q.E.D. 

Finally from this result, one can prove fET]) . For that one realizes (sim- 
ilarly to the necessity proof of theorem 15. ip that a trace and self-adjoint 
preserving generator can be written as 



F^pFl-^{FlFj^P} 



Cp = -i[H,p\ + ^ ttjk 
j,k=i 

where H = H"^ and the coefficients ajk form a Hermitian matrix. 8^- is 
completely positive (not just positive) if and only if £®1 (in the place of C) 
satisfies the conditions ( ISBj) . (|Fr|) and (155]) . It is easy to see [HI] this implies 
that ttjk is positive semidefinite, and then one obtains flSTj) by diagonalization. 
The matrix a^k is sometimes referred to as Kossakowski matrix [28|l83]. 
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5.4 Steady states of homogeneous Markov processes 

In this section we give a brief introduction to the steady state properties of 
homogeneous Markov processes, i.e. completely positive semigroups. This is 
a complicated topic and still an active area of research, so we only expose 
here some simple but important results. 

Definition 5.3. A semigroup Er is relaxing if there exist a unique (steady) 
state pss such that Er (pss) = Pss for all r and 

lim E^ (p) = Pss, 

r— ^oo 

for every initial state p. 

Of course if one state pss is a steady state of a semigroup E^, it is an 
eigenoperator of the operator E^ with eigenvalue 1, 

Er (Pss) ^ Pss Pss; 

and, by differentiating, it is an eigenoperator of the generator with zero 
eigenvalue, 

/:(Pss) = 0. (63) 

If this equation admits more than one solution then the semigroup is ob- 
viously not relaxing, and the final state will depend on the initial state, or 
might even oscillate inside of the set of convex combinations of solutions. 
However, the uniqueness of the solution of equation ([H3D should not presum- 
ably be a guarantee for the semigroup to be relaxing. First of all let us point 
out the next result. 

Theorem 5.3. A completely positive semigroup Er = e^'^ on a finite dimen- 
sional Banach space 53 has always at least one steady state, this is, there 
exists at least a state pss which is a solution of [6l 



Proof. We prove this by construction. Take any initial state p(0); we calcu- 
late the average state of its evolution by (this is sometimes called "ergodic 
average" ) 

p = lim i / e^-'p(0)rfr'. (64) 
It is simple to show that the integral converges just by taking norms 



lim 

T^oo 



\ f e^^pmr < hm ^ T He^^l llp(0) l^r < 1. 
T Jo T-^-^ T Jq 
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Because e'~'^ is completely positive, it is obvious that p is a state. Moreover, 
p is clearly a steady state 

e^V = -p. 

Q.E.D. 

The easiest way to show the steady state properties of a finite dimensional 
semigroup is by writing the generator and the exponential in matrix form. 
In order to do that, let us consider again the Hilbert- Schmidt inner product 
on 03 and a set of mutually orthonormal basis {Fj,j = 1, . . . , A^^}, that is 
(Fj,Ffc)Hs = Tr(Fj^Ffc) = 6jk- Then the matrix elements of the generator £ 
in that basis will be 

Cjk = {Fj,CF,)us = TT[FjC{Fk)]. 

To calculate the exponential of Cjk one uses the Jordan block form (in general 
Cjk is not diagonalizable). That is, we look for a basis such that Cjk can be 
written as a block-diagonal matrix, S^^^S = C^^^ © C^"^^ © ... © 0^''\ of M 
blocks where S is the matrix for the change of basis, and each block has the 
form 

' Xi 1 ... " 

\i 1 ••. ; 
; ••. ••. ••. . 
; ■•. ■•. 1 

Xi _ 

Here Xi denotes some eigenvalue of C and the size of the block coincides 
with its algebraic multiplicity ki. It is a well- know result (see for example 
references [ISIIHI]) that the exponential of C can be explicitly computed as 

e^^ = Se^'^'^ © e^'^'^ © . . . © e^'"'^S-i 

= S (e^i'e^^") © (e^2"e^2") © ... © (e^^^'e^^^") S-\ (65) 

where the matrices Nj = jC^^'^ — Xil are nilpotent matrices, (Nj)^' = 0, and 
thus their exponentials are easy to deal with. If /I is a diagonalizable matrix, 
Nj = for every eigenvalue. 

Now we move on to the principal result of this section. 

Theorem 5.4. A completely positive semigroup is relaxing if and only if 
the zero eigenvalue of C is non-degenerate and the rest of eigenvalues have 
negative real part. 
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Proof. For the "only if" part. The condition of non-degenerate zero eigen- 
value of C is trivial, otherwise the semigroup has more than one steady state 
and is not relaxing by definition. But even if this is true, it is necessary 
that the rest of eigenvalues have negative real part. Indeed, consider some 
eigenvalue with purely imaginary part C{a) = i(f)a, where a is the associated 
eigenvector. By using the fact that e^'^ is trace preserving 



Tr [e^^((T)| = Ilia) = e'^^^Tiia) ^ Tr(a) = 0. 



So we can split a as a difference of positive operators cr = cr+ — cr_, given by 
its spectral decomposition as in the proof of theorem 14.41 Since Tr((T) = 0, 
Tr((j+) = Tr(cT_) = c, and we can write a/c = p+ — p_, where p± are now 
quantum states. The triangle inequality with the steady state pss reads 



\e^-{a/c) 



\e 



Pss) 
- Pss 



Thus if the semigroup is relaxing we obtain lim^ 
is impossible as e^'^{a) = e^'^'^a and 



(a/c) 



Ps 



0, but this 



a/c) 



[a/c) 



\a/c\\ > 0, Vr. 



Conversely, let us denote Ai = the non-degenerate zero eigenvalue of 
C, since the remaining ones have negative real part, according to (165|) in the 
limit T ^ oo the matrix form of the semigroup will be 



lim e^^ = lim S (e^^") © (e^^'e^^") 



)s- 



S(l ©0 



0)S-^ = SBiS\ 



(66) 



where 



Di 



/I 








\o ■■■ y 

coincides with the projector on the eigenspace associated with the eigenvalue 
1. From here it is quite simple to see that the application of £oq on any 
initial state p written as a vector in the basis of {Fj,j = 1, . . . , A^^} will give 
us just the steady state. More explicitly, if pss is the steady state, its vector 
representation will be 



Vss — [(-^1, Pss)hS, • • • , (-^Af2, Pss)hs]* 
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and the condition e^'^p^s = Pss is of course expressed in this basis as 

e^^Vss = Vss- (67) 

On the other hand, as we aheady have said, Di is the projector on the 
eigenspace associated with the eigenvalue 1, which is expanded by Vgg, so we 
can write Di = Vss(vss, ■) or Di = |vss)(vss| in Dirac notation. 

As a particular case (r — oo) of (l67l) . we have ^ooVgs = Vgs, by making a 
little bit of algebra with this condition 

£ooVss = SDiS" Vss = Vss SVss(Vss, S^Vgs) = Vss 

^ SVss = (v..,S-iv..) ^ss, (68) 

therefore Vgg is eigenvector of S with eigenvalue -, — — r. 

l^VssiO Vss j 

Finally consider any initial state p written as a vector Vp in the orthonor- 
mal basis {Fj,j = 1, . . . , A^^}. We have 



lim e^^Vp = £ooVp = SDiS Vp = Svss(vss, S 

T—>-00 

thus using fl68|) . 



fv««,S Vp)__ (Vss,S Vp) 



^ooVp = j -Vss = C'(p)Vss, C{p) = 



(Vss,S IVss) ' (Vss,S V, 

and in the notation of states and operators this equation reads 

lim e'^^p = C(p)pss, for all p. 



T— >00 



One can compute C(p) explicitly and find that C(p) = 1 for every p, however 
this follows from the condition that e^'^ is a trace preserving semigroup. 

Q.E.D. 

This theorem can be taken as a weaker version of the result by Schirmer 
and Wang [85], which assert that if the steady state is unique then the semi- 
group is relaxing (there are not purely imaginary eigenvalues). A particular 
simple case is if the unique steady state is pure pss = \ip){ip\. Then since a 
pure state cannot be expressed as a nontrivial convex combination of other 
states, the equation (l64l) implies that the semigroup is relaxing. This was 
also stated for an arbitrary UDM in theorem 8 of reference [86] . 

On the other hand, it is desirable to have some condition on the struc- 
ture of the generator which assures that the semigroup is relaxing as well 
as a characterization of the steady states; there are several results on this 
topic (see [61 I851I871499] and references therein). Here we will only present an 
important result given by Spohn in [89] . 
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Theorem 5.5 (Spohn). Consider a completely positive semigroup, £r = e'''^ , 
in *B with generator 

Cp = -i[H,p\ + Y,lk \vtpVl - \{VlV,,p} 

for some set of indexes I. Provided that the set {Vk, k & 1} is self-adjoint 
(this is, the adjoint of every element of the set is inside of the set) and the 
only operators commuting with all of them are proportional to the identity, 
which is expressed as {Vk, k G /}' = cl, the semigroup £r is relaxing. 

Proof. Since {Vk,k G /} is self-adjoint we can consider a self-adjoint or- 



thonormal basis with respect to the Hilbert-Schmidt product {Fj,j 



such that, as in the proof of theorem 15. we chose Fi^2 = and the 

remaining elements are traceless. Expanding 



Vfe = > VkmFm H 7=^1 

m=l * 

we rewrite the generator as 

Ar2-1 

Cp = -i[H + H,p\+ J2 

m,n=l 

where H is the self-adjoint operator 

Af2-1 



(69) 



FmPF-n {FfiF^n, 



(70) 



Ik 



o- /Ar {^'k,N2^km - K,N^^'km) F„ 

kel * m=l 



and the Hermitian matrix is a^n = J2kei^kVkmVkn- Now note that we can 
trivially chose amn to be strictly positive; for that it is enough to diagonalize 
again a^n and ignore the members in the sum of fITOl) which goes with zero 
eigenvalue as they do not contribute, and write again the less dimensional 
matrix amn in terms of the numbers of elements of the basis strictly necessary, 
say p < N'^ — 1, 

p 

Cp = -i[H + H,p]+ ^ amn 



m,n=l 



1 



FmpFn 2 {FnFmy p} 



Next, consider some a > to be smaller than the smallest eigenvalue of the 
strictly positive amn, and split the generator in two pieces £ = £i + £2 with 



£2 



a 



E 



m=l 



FmpFm i^{FmFm, p} 
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It is clear that £2 is the generator of a completely positive semigroup, and 
the same for Ci since amn~(^^mn is still a positive semidefinite matrix. Next, 
for some o" G 53 we compute 



(cr,£2(a))HS = TrKsM] =«^Tr[(aF„y-a2i^2 



m=l 

P P 



2 ^ Tr ([F„„ a]) = -- Tr a]t[F„, a]) < 0, 

m=l m=l 

SO £2 has only real eigenvalues smaller or equal to zero. Moreover for the 
eigenvalue we have Tr[o"£2(c")] = which requires = for m = 

l,...,p. But, because of the equation (jSHD, it implies that [Vfc,cr] = for 
k E I and therefore a = cl by presupposition. 

Now, let C, Ci and £2 denote the restriction of £, £1 and £2 to the sub- 
space generated by {Fj.j = 1, . . . , A^^ — 1}, this is the subspace of traceless 
self-adjoint operators. Then since £1 is the generator of a contracting semi- 
group all the eigenvalues of £1 have real part smaller or equal to zero. On the 
other hand from the foregoing we conclude that the spectrum of £2 is com- 
posed only of strictly negative real numbers, this implies that £ = £1 + £2 
has only eigenvalues with strictly negative real part. Indeed, writing £, £1 
and £2 in its matrix representation, C\ and we can split them in their 
Hermitian and anti-Hermitian parts, e.g. 



where Ca = + /2 and Cb = [c - c)^ jli are Hermitian matri- 
ces. Similarly for C\ and £2- It is quite evident that the real part of the 
spectrum of £ is strictly negative if and only if Ca is a negative definite 
matrix, ^v^-, Ca^o^ < for any a G Span{Fj, j = 1, . . . , A^^ — 1}. But, since 

Ca = Cai + Ca2 where Cai = (^Ci + c[^ /2 and £^12 = (^£2 + £2) /2, we 



obtain 



(\^, Ca^J^ = (v^, £AlVa) + (v^, £A2Va) < 0, 



as (^(j, Cai^o^ < and ^v^-, £yi2Vo- j < 0. Thus all the eigenvalues of £ 
have strictly negative real part. 

Finally since the eigenvalues of £ are just the ones of £ plus one (this is 
obvious because the dimension of the subspace of traceless operators is just 
one less than the whole state dimension iV^), and this extra eigenvalue has 
to be equal to because theorem 15.31 asserts that there exists always at least 
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one steady state. Therefore L fulfills the conditions of theorem 15.41 and the 
semigroup e'-''^ is relaxing. 

Q.E.D. 

For checking if some semigroup satisfies the conditions of this theorem 
it is sometimes useful to explore if the operators {V^, A; G /} correspond to 
some irreducible representation of a group. Then we may invoke the Schur's 
lemma [100] . 

6 Microscopic description: Markovian case 

As we have already pointed out, the dynamics of an open quantum system is 
notably different from that of a closed system; both from practical and fun- 
damental points of view. Essential questions as the existence of universally 
valid dynamics or the continuity in time of these, become non-trivial matters 
for the evolution of open quantum systems. 

Generally it would be necessary to describe the evolution in the extended 
closed system, and trace out the state at the end of the evolution. Unfortu- 
nately, many times, one of the systems is out of our control, and/or it has 
an infinite number of degrees of freedom, such that the task of calculating 
the associated Kraus operators can be difficult and impractical. That is the 
case for example of a finite A^- dimensional system coupled with an infinite 
dimensional one; according to Eq. (135|) . we would obtain, in principle, infi- 
nite terms in the sum. However it can be easily shown that any completely 
positive map on a A^- dimensional system can be written with at most A^^ 
Kraus operators (of course; since the dimension of 03* is A^^, so the use of 
the infinite terms is unnecessary. 

For that reason other techniques apart from the "brute force" strategy 
will become very useful and they will be the subject of this section. 

FirstLY, we analyze some conditions under for a UDM to be approxi- 
mately be approximately represented by a quantum Markov process. We 
will start from a microscopic Hamiltonian and apply some approximations 
leading to UDMs which fulfill the divisibility property (HTjl . Note that, in 
principle, as we have already said, the reduced evolution of a quantum system 
is never exactly Markovian. This is because, as discussed in section 1^751 there 
is always a privileged time, say where the both systems A and B start 
to interact being then the global state a product p = Pa® Pb and defining 
a UDM for pA- Since generally the global state can no longer be factorized, 
and (generally) the dynamical map from there is not a UDM, the divisibility 
property (HTj) is violated. 
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6.1 Nakajima-Zwanzig equation 

Let assume that A is our open system while B plays the role of the environ- 
ment; the Hamiltonian of the whole system is given by 

H = Ha + Hb + V, (71) 

where Ha and Hb are the Hamiltonians which govern the local dynamics of 
A and B respectively, and V stands for the interaction between both systems. 

We will follow the approach initroduced by Nakajima [lUlj and Zwanzig 
|1U2] using projection operators, which is also explained in [SlISniEI] among 
others. In this method we define in the combined Hilbert space "system plus 
environment" "H = 'Ha^'Hb, two orthogonal projection operators, V and 
g, = V^ Q"^ = Q and PQ = QV = 0, given by 

Vp = Ttb{p)®Pb, (72) 
Qp = il-V)p, (73) 

where p G "H is the whole state of the system and environment and pb &'Hb 
is a fixed state of the environment. In fact, we choose pb to be the physical 
initial state of the environment. Strictly speaking this is not necessary but 
the method becomes much more complicated if one chooses another state. 
In addition, we assume that the system B is initially in thermal equilibrium, 

tills IS 

PB = Pth = e-^^-[Tr(e"^^-)]-\ 

where /3 = 1/T. Note that Vp give all necessary information about the 
reduced system state pA', so the knowledge of the dynamics of Vp implies 
that one knows the time evolution of the reduced system. 

If we start from the von Neumann equation (1191) for the whole system 



J^ = ^,IHA + HB + V,p{t)], (74) 
and take projection operators we get 

Ipp(t) = -^r[HA + HB + V,p{t)], (75) 

^Qp(t) = -iQ[HA + HB + V,p{t)]. (76) 



As usual in perturbation theory we shall work in interaction picture with 
respect to the free Hamiltonians 

p{t) = e^(^^+^s)V(t)e-'(^-^+^^)*. 
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Since 

TTB[p{t)] = TiB [e'(f^A+Hs)tp^^^^-iiHA+Hs)t^ 

the projection operators are preserved under this operation, and equations 
fl75|) and fl76l) can be written in interaction picture just hke 

^Vp{t) = -^V[V{t)rp{t% (77) 
j^Qp{t) = -zQ[V{t)rp{t)]. (78) 

Let us use the notation V{t)- = —i[V{t), ■]; then by introducing the iden- 
tity 1 = V + Q between V{t) and p{t) the equations ( 177|) and (1751) may be 
rewritten as 

±Vm = VV{t)Vp{t) + VV{t)Qp{t), (79) 
^Qp(t) = QV{t)Vp{t) + QV{t)Qp{t). (80) 
The formal integration of the first of these equations gives 

Vp(t)=Vp{<d)+ [ dsVV{s)Vp{s)+ [ dsVV{s)Qp{s), (81) 
Jo Jo 

where we have set as the origin of time without losing generality. On the 
other hand, the solution of the second equation can be written formally as 

Qpit) = G{t,0)QpiO)+ [ dsg{t,s)QV[s)Vp{s). (82) 

Jo 

This is nothing but the operational version of variation of parameters formula 
for ordinary differential equations (see for example |1031I104] ). where the 
solution to the homogeneous equation, 

jfip{t) = Qvmpit). 

is given by the propagator 
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Inserting ([12]) in ([HI]) yields 

rp{t) = Vp{0)+ [ dsVV{s)Vp{s) + ! dsVV{s)g{s,^)Qp{Q) 
Jo Jo 



ds / duPV{s)g{s,u)QV{u)Vp{u). 
Jo 



(83) 



This is the integrated version of the so-called generalized Nakajima-Zwanzig 
equation 



d_ 
di 



vp{t) = vv{t)rp{t) + vv{t)g{t,o)Qp{o) 



duVV{t)g{t,u)QV{u)Vp{u) 



Next, two assumptions are usually made. First one is that VV(t)V = 0, 
this means, 



VV{t)Vp = -i TiB[V{t),pA ® Pb] ® pB 



TrB ( V{t)pB ] , Pa 



® = 0, 



for every p and then every p^, which implies Tr^ {v{t)pB^ = 0. If this is not 

fulfilled (which is not the case in the most practical situations), provided that 
Pb commutes with Hb one can always define a new interaction Hamiltonian 
with a shifted origin of the energy for A [TOlIM]. The change 

V' = V -TtbIVpb]®!, and H'^ = Ha + Ttb[Vpb] ® 1, (85) 
makes the total Hamiltonian to be the same and 



Tr 



B 



V'it)p 



B 



Tr 



B 



— e 



TiB[VpB]e~'"'^'TiB[e'''^'pBe 



Hut. 



-iHgt] 



e'''^'TTB[VpB]e- 



JH', t I 



e-^'TiB[VpB]e-'''^' = 0, 



as we wished. 

The second assumption consists in accepting that initially p(0) = pa(0) Cg> 
Pb(0), which is the necessary assumption to get a UDM. Of course some 
skepticism may arise thinking that if the system B is out of our control, there 
is not guaranty that this assumption is fulfilled. Nonetheless the control on 
the system A is enough to assure this condition, for that it will be enough to 
prepare a pure state in A, for instance by making a projective measurement 
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as argued in |1U5] . to assure, by theorem 14. II that at the initial time the global 
system is in a product state. 

These two assumptions make to vanish the second and the third term in 
equation (183|) . and the integro-differential equation yields 



If this kernel is homogeneous Kit, u) = }C(t — u), which is the case for station- 
ary states pb (as exemplified in next section), this equation can be formally 
solved by a Laplace transformation but this usually turns out into an 
intractable problem in practice. 

In order to transform this integro-differential equation in an Markovian 
master equation one would wish that the integral on the right hand side of 
fISB]) turns into a generator of the form of fH^ . For that aim, the first intuitive 
guess is to try that /C(t, u) behaves as a delta function with respect to p{u). 
For that to be true the typical variation time ta of p{u) (which is only due 
to the interaction with B because we have removed the rest of the dynamics 
by taking integration picture) has to be much larger than some time tb, 
which characterizes the speed at which the kernel /C(t, u) is decreasing when 
|t — m| ^ 1. Of course this kind of approximations intrinsically involve some 
assumptions on the size of B and the "strength" of the Hamiltonians, and 
there is in principle no guaranty that the resulting differential equation has 
the form of fH6|) . However, several works |106H110] studied carefully two 
limiting procedures where ta/tb — ?■ oo so that this turns out to be the case. 
These two limits are 

• The weak coupling limit. Let us rewrite V — > aV, where a accounts 
for the strength of the interaction. Since for a small the variation of 
p{u) is going to be slow, in the limit a — > one would get ta — )• oo, 
and so ta/tb oo for every fixed tb- Of course for a = there 
is not interaction, and so we will need to rescale the time parameter 
appropriately before taking the limit. This is discussed in more detail 
in next section 16.21 

• The singular coupling limit. This corresponds to somehow the oppo- 
site case; here we get — )■ by making /C(t, u) to approach a delta 
function. We will discuss this limit in the forthcoming section 16.31 




(86) 



where 



/C(t, u) 



VV{t)g{t,u)QV{u). 
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We should stress that one can construct Markovian models ad hoc which 
try to describe a system from a microscopic picture without involving these 
procedures; see for example 1111(1112] . In this sense, these two cases are 
sufficient conditions to get a Markovian evolution, but not necessary. 



6.2 Weak coupling limit 

As we said above, if we assume the weak coupling limit, the change of p{u) 
is negligible within the typical time tb where the kernel /C(t, u) is varying. 
To perform this limit, consider again Eq. f l83|) under the assumptions of 
VV{t)V = and p(0) = Pa{0) ® ps(0): 



Vp{t)=Vp{0) + a^ [ ds [ duVV{s)g{s,u)QV{u)Vp{u) 
Jo Jo 



i7) 



where we have redefined V — aV (so V — )■ aV). Since in the interaction 
picture the whole state satisfies a von Neumann equation of the form 



d 
dt 



p{t) = -ia[V{t),p{t)l 



we can connect the whole state at times s and u (note u < s) hj the unitary 
solution 

p{s) = U{s,u)p{u)U'^ {s,u) = U{s,u)p{u) 

^ p{u) = W{s, u)p{u)U(s, u) = U{u, s)p{s). 

By introducing this operator in the last term of fISTI) we find 

Vp{t)=Vp{0)+a^ [ ds f duVV{s)g{s,u)QViu)VU{u,s)p{s). 
Jo Jo 

Now the term inside the double integral admits an expansion in powers of a 
given by the time-ordered series of 

g(s,u) = Te'^n'it'QVit') ^^^^ 

U{u,s) = [W(.,n)l^= fre°/^^'^*'^(*'f = r'^e"/^^'^*'^'(*'), (89) 



where the adjoint is taken over every element of the expansion and T* denotes 
the antichronological time-ordering operator. This is defined similarly as T 
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in f[T^ but ordering in the opposite sense, this is, the operators evaluated at 
later times appear first in the product. Thus, at lowest order we get 

Vpit) = Vp{0) + a^ [ ds ! duVV{s)QV{u)Vp{s) + 0{a^) 



^0 

t 



= Vp{Q) + a^ f ds f duVV{s)V{u)Vp{s) + 0{a^). 
Jo Jo 

where we have used again that VV{t)V = 0. After a change of variable 
u — )■ s—u, writing explicitly the projector V and the operators V, we conclude 
that 

pA{t) = Pa(0) -a^ [ ds f du TiB \v{s), [V{s - u)rpA{s) ® Pb\ \ + 
Jo Jo ^ 

(90) 

6.2.1 Decomposition of V 

In order to continue with the derivation, we write the interaction Hamiltonian 
as 

V = Y,Ak®Bk, (91) 

k 

with A\ = Ak and 5^ = Bk- Note that since V'^ = V, one can always chose 
Ak and B^. to be self-adjoint. Indeed, assume that ^ is a sum of products of 
any general operators and Y^, it will be written then like 

n 

y = ^ ® n + Xfc ® Y^. (92) 

k=l 

Any operator can be decomposed as 

X, = Xi^^ + ^X!'\ 

where 

(,) _ {Xj + X,) (,) _ ^iXl-X,) 
- 2 ' ' ~ 2 
are self-adjoint. This, and the same kind of decomposition for Yf^ = Y^'^^ + 
lY^^ inserted in equation (192|) gives 

n 

V = 5:(4">-af')«(yi"' + ,yf' 

n 
k=l 
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fc=i 



so one can make the choice 

Ak = 2X\^\ Bk = 2Ylf\ for A; = 1, . . . , n 

Ak = 2Xf2„, 5fc = 2Fi!)„, for A: = n + l,...,2n 

to arrive at ( 19T]) . 

Now we focus on the operators A^, let us assume for the moment that 
the spectrum of Ha is discrete, and let \ips) be an eigenstate associate with 
the eigenvalue e. We define 

Ak{uj)= \^IJe){^IJe\Ak\i!e'){llJe'l 
e' —£=iij 

where the sum is over every e' and e such that their difference is u. The 
operator so defined is an eigenoperator of C-^ = i[HA_, ■] with eigenvalue —iu; 
indeed, 

C^Akiu) =i[HA,Akiuj)] = t {e-e')\ij,){^lj,\Ak\^IJe'){A'\ = -^(^M^)- 

e'—e=u] 

As a result, 

C-^Aliu) = tuAliu), 

with 

and this holds for every k. In addition, it is straightforward to verify that 

[HA,Al{u)Ae{io)] = 0. (93) 

On the other hand, to write A^^u) and in the interaction picture 

with respect to Ha is quite easy: 

e^^'Akiu) = e'"^'Ak{u)e-'^^' = e-''^'Ak{u;), 
e^^'Aliu) = e'''^'A\{ujy-'''^' = e'^'A\{uj). 

Moreover, note that 

^ AfeH = Y \'^e){i'e\Ak\^e')hl^e'\ = ^fc, 
ui e' ,£ 

and similarly that 
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Thus, by introducing these two decompositions for Aj^ in (1911) and taking the 
interaction picture, one gets 

V{t) = J2 e-'^^'Akiuj) ® Bk{t) = J2 e'"*4(^) ® ^IW- (94) 

ui,k u),k 



Now, coming back to the equation fl90|) the use of the above first equahty 
with V{s — u) and the second one with V{s) yields 

PA{t) = pA{0) + rrfs5^5^e^(-'--)T^,M[A,Mp^(.),4(a;')] 

uj,uj' k,l 

+ e'(— ')T-(a;)[A,(a;'),PA(s)4(^)] + (95) 
Here we have introduced the quantities 



due*'^"Tr 



Bk{s)Bi{s - u)pb 
Bk{u)BiPB 



(96) 



where the last step is justified because pb commutes with e^^'^^. Next, we 
take a rescaled time r = a'^t and a = a'^s, obtaining 

Pa(^) = pAir/a"^) 

= pa{0)+ rda5^5^e^(-'--W"^r^/"'(a;)[yl,(a;)p^(a/a^),4(^')] 



+ e 



i(u]-uj')a/a'^-p^/a'^ 



Ik 



{uj)[A,{uj')rPA{ala^)Ai{uj)] + 0{a). 



From this equation we see that in the limit a — )■ (keeping finite r and a) 
the time scale of change of pAij/a^) is just r. This is merely because the 
dependency with of pa[jIo?^ enters only in the free evolution which we 
remove by taking the interaction picture. The remaining time scale due to 
the interaction with the environment is r, as expressed in the above equation 
(the integral runs from to r), so we write 



limp^(t) = lime^-^^"/°VA(r/a^)e-^-f^^"/°' = pAij) 

«— >-0 a— >-0 
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and thus 



Pa{t) = p^(0) + lim rrfa5^5^e^(-'--W°^r^/"'M[A,(c.)p^(a/«^),4(c.')] 

^0 . . . .1 u c 



-,,2 



Before going further in the derivation, a mathematical result which will 
be particularly useful from now on is the following. 

Proposition 6.1 (Riemann-Lebesgue lemma). Let f(t) be integrable in [a,b] 
then 



lim / e^^*/(t)rft = 0. 

J a 



Proof. There are several methods to prove this, here we define g{t) = f{t) [6{t- 
a) — 9{t — b)] and so 



h poo 

ixt rfj.\ 7i / ixt , 



e'^'f{t)dt= / e'^'g{t)dt = gix) 



oo 



where g{x) denotes the Fourier transform of g{t). Now it is evident that g(t) 
is square integrable, then as the Fourier transform preserve norms [IS] (this 
is sometimes referred to as Parseval's theorem): 

/oo 
\gix)\''dx, 
'OO 

which implies that g{x) is also square integrable and then lima-_^oo .^(2;) = 0. 

Q.E.D. 

Because of this result (applied on the Banach space to be concise) 
in the limit a — ?■ the oscillatory factors in ( 197|) with u' ^ u are going 
to vanish. This so-called secular approximation is done in the same spirit 
that well known rotating wave approximation in the context of quantum 
optics [SOIEIIESIEHIEI]. We obtain 



^0 k,i 

= pa{0)+ [ daC[pA{a)]. (9^ 
^0 



60 



Here the coefficients are given by the one-sided Fourier transform 



Bk{u)BipB 



6.2.2 Correlation functions 

In order to make convergent the above integrals, the environmental correla- 
tion functions Tr Bk{u)BipB have to decrease as u is increasing. However, 
if we use the corresponding eigenoperator decomposition of Bk for the Hamil- 
tonian Hb-, we find 



Tr 



Bu{^i)B,pB] =^e-^""Tr[5fc(u;)5,pB] 



which shows that the correlation functions are periodic in u and their inte- 
gral extended to infinity will diverge. The only possible loophole to break 
the periodicity is the assumption that the system B has infinite degrees of 
freedom, in such a way that Hb has a continuous spectrum. Then every- 
thing is essentially the same but the sums are substituted by integrals, i.e. 
the decomposition in eigenoperators will be 



Bh 



dujBkioj) 



(99) 



where a is the maximum eigenfrequency (it can be infinite). So the correlation 
functions become 



Tr 



Bk{u)BipB 



due-'^''Ti[Bj,{u)B,pB] 



and of course now they are not periodic as 



lim Tr 

because of the proposition (16. 



Bk{u)BipB 



6.2.3 Davies' theorem 

We have just seen that a necessary condition for the weak limit to exist is that 
the environmental system B has infinite degrees of freedom. However this is 
not a sufficient condition because the decreasing behavior of the correlation 
functions may be not fast enough for the one-sided Fourier transform to exist. 
A sufficient condition was pointed out by Davies. 
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Theorem 6.1 (Davies). // there exist some e > such that 



dt 



Tr 



Bk{t)BepB {l-ty<oo, 



(100) 



then the weak coupling is strictly well defined (for a bounded interaction V) 
and 

lim ^ ||4,o) [pa(0)] - e^" [pa(0)] ||i = 

for all Pa(0) and uniformly in every finite interval [0, tq], where £{Tfi) denotes 
the exact reduced dynamics in the interaction picture. 

Proof. To prove the condition fllOOl) is beyond the scope of this work, it was 
developed in |106ill07] . Under f llOOp the convergency to the above result is 
clear from equation (j98p which is the integrated version of 

dpAir) 



dr 



PA{T} 



[PAm ■ 



Q.E.D. 



On one hand, note that the method works because there are no eigen- 
frequencies, w, arbitrary close to one another, as a result of the discreteness 
of the spectrum of H^, otherwise the secular approximation is jeopardized. 
For the infinite dimensional case where this is not satisfied, recently D. Taj 
and F. Rossi |113] have proposed another approximation method which sub- 
stitutes the secular approximation and also leads to a completely positive 
semigroup. 

On the other hand, an important consequence of the secular approxi- 
mation is that it guaranties that the generator C has the correct form of 
a Markovian homogenous master equation (theorem 15. ip . To see this, let 
us decompose the matrices r^(co') in sum of Hermitian and anti-Hermitian 
parts 

(101) 



where the coefficients 



and 



Mc^) = ^[rSM-r,TM], 



due-''^''Ti 



Bk{-u)Bepth 

(102) 
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form Hermitian matrices. In terms of these quantities, the generator L can 
be written as 

t\pA{T)\= - z[i/LS, PA(r)] 

+ IZIZ^fc^M ^t{^ypA{r)A\{uj)-\-{A\{uj)A,{uj)rPA{r)') 



UJ k,i 

where the Hamihonian part is given by 



^LS = 5Z 5^ SU^)Al{u)Mto). (103) 

In order to prove that C has the form of a generator of a Markovian homoge- 
nous evolution, it is required just to show that the matrix •jkei^) is positive 
semidefinite for all u. This is a consequence of the Bochner's theorem [^H5]. 
which asserts that the Fourier transform of a function of "positive type" is a 
positive quantity. Note that a function f{t) is of positive type if for any 
and tn the matrix constructed as fmn = fitm — tn) is positive semidefinite. 
For any vector v we have 

/oo 
(iwe"^" vlTi Bk{u)Bipth ve 
K,t k,£ 

/oo 
k,£ 

/oo 
-oo 

where C = Xlfc BkVk- 

Proposition 6.2. The function f{t) = Tr [e^^BtQl^-iHBt^jp^^-^ "^^g^. 
tive type" and then the matrix is positive semidefinite (v,7v) > 0. 

Proof. The positivity of (v, 7v) > follows from the Bochner's theorem if 
f{t) is of "positive type" because it is its Fourier transform. For the proof of 
the first statement, let us take the trace in the eigenbasis of Hb (remember 
that the spectrum is taken to be continuous): 

fit) = j rf£((p,|e^^^*Cte-^^«*Cpth|v'e) = j &e-Ve(v'.|Cte"^^^*C|^,), 



here we have used that 
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with Pe > since pth is a positive operator. By introducing the identity 

1 = J de'\(pe'){(Ps'\, 

fit) = J dede'e'''p,{if,\C^e~'"-'\^,>){if,>\C\^,) 
dede'e'^^-^'^%\{^AC\Ve)?. 



Finally we take another arbitrary vector w and the inner product 

(W, fw) = W*^f{tin - tn)Wn 

m,n 



dede' 



E 



because the integrand is positive on the whole domain of integration. 

Q.E.D. 

As a conclusion of these results, if for a small, but finite a, we substitute 
the exact reduced dynamics S(^t,o) by the Markovian one e^"^ = e° the error 
which we make is bounded and tends to zero as a decrease provided that the 
assumptions of the Davies' theorem are fulfilled. 

Of course one can write a'^C and immediately go back to Schrodinger 
picture by using the unitary operator e*^-** to obtain 



dpAjt) 
dt 



(104) 



U) kl 



A,{uj)pA{t)Al{uj) - ]^{Al{uj)A,{uj),pA{t)] 



Note that because of (l93i) . the Hamiltonian i^LS commutes with Ha-, so it 
just produces a shift in the energy levels of the system A. Thus this is 
sometimes called Lamb shift Hamiltonian. However, note that i^LS is not 
only influenced by an environment in the vacuum state (zero temperature) 
but also by thermal fluctuations. 
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6.2.4 Explicit expressions for and 5*, 



kt 



We can give explicit expressions for 7^^ and Sm by using the decomposition 
of the operators Bk in eigenoperators. For the decay rates 

/oo _ 
rfwe^-^Tr Bk{u)B,p,^ 
-00 

/a poo 
dto' / due'^''-'^>TT[Bkico')B,p,^] 
-a J —00 

= 2n f du'6iu-u')Tr[Bkiuj')BePth] 

J —a 

= 27iTT[Bk{u)B,pt^], (105) 

provided that u G (—a, a). So •jkei,^) is just proportional to the correlation 
function of the eigenoperator Bkioj) with the same frequency oj, unless this 
frequency is outside of the spectrum of [Hb-, ■], then 7fc^(a;) will be zero. 

To get an explicit expression for the shifts Ski. is a little bit more compli- 
cated. From fllOll) we have 

1 

Ski{uj) = -ir^(w) + -7A,.f(w). 

The substitution of the decomposition in eigenoperators in r^^(co') and the 
previous result fll05p for 7fc,^(w) yields 

/a POO 
dcu' / due'^'^-'^^Ti [Bkiuj')BepB] + vrzTr [5fc(a;)fi,pth] • 
-a Jo 

(106) 

Note that the integral due^^'^~'^ ^" in not well-defined. Similarly it hap- 
pens with dwe*^'^"'^'^", but we know this integral is proportional to the 
delta function. The mathematical theory which provides a meaningful in- 
terpretation to these integrals is the theory of distributions or generalized 
functions. Here we will not dwell on it, references ^011^ provide introduc- 
tions and further results can be found in [114(1115] . For us it is enough to 
show that these integrals can have a meaning if we define them as a limit of 
ordinary functions. For instance if we define 

poo roo pO 

I due'^'^-'^> = hm / due'^'^~'^>-''' + lim / due*(^^'^')"+™, 

^^0+ Jo e^O- J_^ 
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by making the elementary integrals, 



due'^'^-'^> = lim ^ : + lim 



e^0+ [U — UJ') + ei e^O- €1 — [UO — LU') 

lim 



e^o (cj — uj') + ei ei — [uj — cu') 
= lim^— -^^ — = 2'jt6{u - u'). 

We identify the well-known family of functions f^{uj—uj') = ^^^24.(^_^/)2] which 
tend to the delta function as e — j- |117] because for any (regular enough) 
integrable function g{u) 

lirn / duj' f\{ijj — uj')g{uj') = g{uj), {a < u < b). 



a 



So, if we do the same procedure but just with the one-sided Fourier transform 
Jg°° due^^^~^ we obtain in equation fll06p 

r 1 

Sm{,uj) = hm / du'- — — :Tr [Bkiuj')BepB] + vrzTr [Bk{uj)Bepth] ■ 

Multiplying the integral by the complex conjugate of the denominator we 
obtain 

j.^ r^^. Tr[i?,(cy)^,p^] ^ r du:'^^Tr[B,{u:')B,p,] 

- £m i j duj \^_J^^^^^ i:i [Bk{uj')BtPB\ 

/a ^ 
^^' (o^-Jr^+e^ Tr [Bk{uj')BtpB\ 

- iTxTi [Bk{uj)BipB] ■ 
Therefore the last term is canceled in the shifts, which read 

/a ^ 

Now we decompose the integration interval {—a, a) in three subintervals 
(—a, u — 6), {u — 6,uj + 6) and {u + S, a) and take the limit 5 — 

^^'oI^TpT^Tr [Bk{uj')BipB] 
+ lim / duj'r^^^^^^T^[Bkioo')BipB] 
+ hm / du \J:^-^:i., TT[B,{u;')B,pB]. 
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The last integral just goes to zero because for any e as small as we want 
the integrand is an integrable function without singularities. Whereas, since 
lime_s>o+ (a;^)2+e2 ~ {lj~lu') ' limiting process in 5 is just the definition of 
the Cauchy principal value of the integral |117] : 



^-^0+ [J^^ {u-uj') J^+s {uj-uj') 

^ r^,iWM^. (107) 

J -a - W') 

In theory of distributions, this result is sometimes referred as Sochozki's 
formulae |115] . 

6.2.5 Typical examples 

It is not worth to dwell on examples of the application of the above results, 
as in the books are plenty of them. Just for sake of comparison, we sketch 
some of the most common cases. 

• Two-level system damped by a bath of harmonic oscillators. The to- 
tal Hamiltonian is 

H = Hsys + -ffbath + V 

= -^(^z + / duala^ + / dujh{uj) {a+a^ + a_al) , 



where is the corresponding Pauli matrix, a± = {a^ ± iay)/2, and 
a^j accounts for the continuous bosonic operators [atj,a^,] = 6{uj — 
u'); where u is the frequency of each mode. The upper limit uj^ax is 
the maximum frequency present in the bath of harmonic oscillators, 
it could be infinity provided that the coupling function h{u) decreases 
fast enough. We have taken the continuous limit from the beginning, 
however it is common to do it in the course of the computation. 

The interaction Hamiltonian may be written in the form V = Ai ^ 
Bi + A2 ^ B2, where it is easy to check that the system operators and 
their eigendecompositions are given by 

Ai = aa: = cr+ + a_, Ai{uo) = cr^, Ai{-UJo) = cr+, 

A2 = (Ty = - cr_), ^2(^0) = -icT-, A2{-uJo) = ia+. 

Similarly for the bath operators, 

Bi= dujh{uj)^^^^^ = / dujBi{uj), 
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with Bi{u) = Bi{-uj) = and 



B2= dujh{uy-^^^^= / dujB^iuj), 

JO J — C^max 

with B2iuj) = zihM^^ B^i-cu) = 

The decay rates are given by 'jk/i^o) = 27r TT[Bk{ujQ)BipB], particularly 
7ii(^o) = 27rTi[Bi{ujo)BipB] 

= 2^*^^°'' j di^jh{u)Ti[a^^X(^^ + al)pB\ 

TT TT 

Here J{ujq) = h'^iuJo) is the so-called spectral density of the bath (which 
accounts for the strength of the coupling per frequency) at frequency 
Uo and n{ujo) is the mean number of bosons in the thermal state pB of 
the bath with frequency cuq; i-e. the expected number of particles of 
the Bose-Einstein statistics 



where T is the temperature of the bath. For the rest of the elements 
one obtains essentially the same 

iM = ^J{uJo)[n{uJo) + 1] 1 
Similarly for 'jkA-^o), 

Let us make the analog for the shifts 5*^^; the first element is given by 

,Ti[B,{uj')B,pb] 



/i^max 
dio 



max 
t^max 



P.V. / du 



{uJo - LO') 

^, TT[Bk{uj')BePB] 
^" (^^0 - t^O 

, Tr[Bki-uj')BePB] 
(^0 + (^') 

n{uj') + 1 n{u') 



^t^max 

P.V. / du 



^max 



^P.V. / du'J{u' 
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_ {uo - u') {uo + cu') _ 



and the rest of the elements are 
5*22(^0) = 5*11 (wo), 

S12M = 5*iM = tP-V. 



duj'J{u') 



n[uj 



n{uj') 



For the shifts with —ujq: 

Sii{-UJo) = S22{-UJo) 

5*12 (-Wo) = S^ii-ujo) 



— P.V. j du'J{u') 



n[uj' 



4 



-P.V. 



duj' J iu)') 



n(u') + 1 



_ {uJo + uj') {uJo - w') 
n{uj') + 1 n{uj') 



Therefore according to fll03p the shift Hamiltonian is 

Jo (^0 - w ) Jo (c^o - ) 



A' + y I a. 



where A' is the integral depending on 77,(0;') which leads to a shift due 
to the presence of the thermal field (Stark- like shift), and A is the 
second integral independent of 77(0;') which lead to a Lamb-like shift 
effect. Of course the identity 1 does not contribute to the dynamics as 
it commutes with any p, so we can forget it. 

For the total equation we immediately obtain 



dt 



Wo + A 



+ A' a^,psys{t) 



+ T[n{coo) + 1] 



r77(wo) 



(T^Psys{t)(r+ - -{a+a-,Psys(t)} 



1 



where F = 2-11 J (uq). 

Harmonic oscillator damped by a bath of harmonic oscillators. This case 
is actually quite similar to the previous one, but the two-level system 
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is now substituted by another harmonic oscillator: 



H 



— Hsys + i/bath + y 



= Uoa'a + 

Jo 



dujal,a^ + 




The bath operators Bi and B2 are the same as above and the system 
operators are now 

Ai = + a, Ai{uo) = a, Ai{—uo) = a\ 

A2 = i{a^ - a), A2{uo) = -ia, A2{-uo) = iaK 

Thus by mimicking the previous steps we obtain 



Note in this case that, on one hand, the Stark-like shift A' does not 
contribute because of the commutation rule of a and a^; and, on the 
other hand, the dimension of the system is infinite and the interaction 
with the environment is mediated by an unbounded operator. How- 
ever as long as domain problems do not arise one expects that this is 
the correct approximation to the dynamics (and actually it is, see for 
instance [36]). 

• Pure dephasing of a two-level system. In this case the interaction com- 
mutes with the system Hamiltonian and so the populations of the initial 
system density matrix remain invariant. 



This is an interesting case to explore because it involves some subtleties, 
but it is exactly solvable (see for example chapter 4 of [5T], |116] . and 
references therein). The interaction Hamiltonian is decomposed asV = 
A ^ B, where A = A{u = 0) = az (o"^ is eigenoperator of [a^, ■] with 
zero eigenvalue!) and 



dt 



= -i[{uJo + A)az,Psys{t)] 

+ r[n(a;o) + 1] apsys{t)a^ - ^{a^a, psys{t)} 

+ Tn{uo) Psys{t)a - l-{aa\ psys{t)} . 



H 



Hsys + -f^bath + y 





for UJ > 0. 
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So some problem arises when we calculate the decay rates because 
B[u} = 0) is not well-defined. The natural solution to this is is to 
understand u = limit in such a way that 

7(0) = 2iT hmTi[B{co)BpB], 

but then there is a further problem. Depending on which side we ap- 
proach to the function Tt[B{u)Bpb] takes the value 27rJ(a;)[n(a;) + 1] 
for cu > or 27rJ(|a;|)n(|a;|) for a; < 0. Thus the limit in princi- 
ple would not be well-defined, except in the high temperature regime 
where ~ "/^(li^l) + 1- However both limits give the same result 

provided that one assumes lim^^o '-'^(li^l) = 0. This is a natural phys- 
ical assumption because a mode with frequency does not have any 
energy and infiuence on the system. Thus if lim^^^o •^(I'^l) = we find 

7(0) = 27r \un J{\u\)n{\u\). 

Actually, since goes to infinity as l/\uj\ when u decreases, if 

J(|ci;|) does not tend to zero linearly in u (this is called Ohmic-type 
spectral density [33]) this constant 7(0) will be either or infinity. This 
fact restricts quite a lot the spectral densities which can be treated 
for the pure dephasing problem in the weak coupling framework with 
significant results (which is not a major problem because this model is 
exactly solvable as we have already pointed out). On the other hand, 
shifts do not occur because i^LS oc CzCTz = 1, and finally the master 
equations is 

+ 7(0)[ (t)az - Psys(t)] • 

6.2.6 Some remarks on the secular approximation and positivity 
preserving requirement 

As we have seen, the weak coupling limit provides a rigorous mathematical 
procedure to derive Markovian master equations. Nonetheless, in solid state 
or chemical physics the interactions are typically stronger than for electro- 
magnetic environments, and it is sometimes preferred to use the perturba- 
tive treatment without performing the secular approximation. However, this 
method jeopardizes the positivity of the dynamics [118] (see also [10511119] ) 
and several tricks have been proposed to heal this drawback. For exam- 
ple, one possibility is taking in consideration just the subset of the whole 
possible states of the system which remain positive in the evolution |12Uj . 



—I 



— az,psys[t) 
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or describing the evolution by the inclusion of a "slippage" operator [121] . 
These proposals can be useful just in some situations; in this regard it seems 
they do not work well for multipartite systems |122[ll23j . 

We have already explained in detail the requirement of the complete pos- 
itivity for a universal dynamics. The option to proceed without it may be 
useful but risky, and actually the risk seems to be quite high as the positivity 
of the density matrix supports the whole consistency of the quantum theory. 
We believe the secular approximation should be used as it is mathematically 
correct in the weak coupling limit. In case of the failure of a completely 
positive semigroup to describe correctly the dynamics of the physical sys- 
tem, it is probably more appropriate the use of non-Markovian methods as 
the "dynamical coarse graining" (see section 17. 3p which preserves complete 
positivity. 



6.2.7 Failure of the assumption of factorized dynamics 

In some derivations of the weak coupling limit it is common to use the fol- 
lowing argument. One makes the formal integration of the von Neumann 
equation in the interaction picture (l88l) 



to give 



p(t) = p(0)-za / ds[V{s),p{s)]. 
Jo 

By iterating this equation twice one gets 

ft, pt P S 

p{t)=p{0)-ia / ds[V{t),p{0)]- ds du[V{s),[V{u), p{u)]]. 







After taking partial trace, under the usual aforementioned assumptions of 
p(0) = pyi(O) ® Pb(0), and Tr^ V(t)pB{0) = 0, this equation is simplified 
to ^ 

pA{t) = Pa{^) - I ds ! dui:iB[V{s),[V{u)rp{u)]]. 
Jo Jo 

Now, one argues that the system B is typically much larger than A and so, 
provided that the interaction is weak, the state of B remains unperturbed by 
the presence of the A, and the whole state can be approximated as p{u) ~ 
Pa{u)® Pb{u). By making this substitution in the above equation one arrives 
at equation fl90l) (modulo some small details which one can consult in most 
of the textbooks of the references), and henceforth proceed to derive the 
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Markovian master equation without any mention to projection operators, 
etc. 

This method can be considered as an effective model in order to obtain 
equation ( 190|) . however it does not make any sense to assume that the physical 
state of the system factorizes for all times. Otherwise how would it be possible 
the interaction between both subsystems?. Moreover, why not making the 
substitution p{u) ~ Pa{u) ® Psiu) directly in the von Neumann equation 
flHHj) instead of in its integrated and iterated version?. 

The validity of the substitution p{u) ~ pa{u) ® Pb{u) is justified because 
of the factorized initial condition and the perturbative approach of the weak 
coupling limit. Of course we pointed out previously that the weak coupling 
approach makes sense only if the coupling is small and the environment has 
infinite degrees of freedom. This fits in the above argument about the size 
of 5, but it can dismissed even by numerical simulations (see [36]) that 
the real total state p{u) is close to pa{u) ® pb{u). It should be therefore 
be considered as an ansatz to arrive at the correct equation rather than a 
physical requirement on the evolution. 



6.2.8 Steady state properties 



In this section we shall analyze the steady state properties of the Marko- 
vian master equation obtained in the weak coupling limit f ll04p . The most 
important result is that the thermal state of the system 



Tr[e(-^A/T)] ' 



with the same temperature T as the bath is a steady state. To see this one 
notes that the bath correlation functions satisfies the so-called Kubo-Martin- 
Schwinger (KMS) condition [124] . namely 



Tr 



PthBk{u)Bi 
1 



1 



■Tr [e'"^^''+''^^Bke-'^^''Bi\ 



Tr [E^e^^^("+^'')5, 



-iHgiu+ip) 



pth] = {B,Bk{u + iP)). (108) 
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By decomposing in eigenoperators B{uj') of [Hb-, ■], we obtain, similarly 
to fllOSp . the following relation between their Fourier transforms 



rfne^'^^Tr Bk{u)B,p,^ 

-oo 
oo 

rfne^'^"Tr B,Bk{u + i(5)p,^ 

-oo 

a POO 

dw'e^"' / due'^'^^'^>TT[BeBk{u')pth] 

-a J —oo 

2iT I du'e^^' 5{uj - uj')i:i[B,Bk{uj')p,^] 

J —a 

27re^"Tr [B,Bk{uj)p,^] = 27re^" jxr \bI{uj)B,p,^ 
27re^-{Tr[5fc(-a;)5,pth]}* 



(109) 



On the other hand, since A{uj') is an eigenoperator of [if^,'], one easily 
obtains the relations 



(110) 
(111) 



Now one can check that = according to equation (11041) . Indeed, it is 
obvious that commutes with the Hamiltonian part, for the remaining one 
we have 



d^ 
dt 



UJ 



Mu)p'}Al{uj) - ^{Aliuj)A,{u),p'^} 



ui kl 



e -^-A,{u;)Al{u;)p'^ - Aliu;)A,{u;)p'^ 
[lki{co)e-^^Al{-co)Ak{-co) - 7,, (a;) 4 (a;) A, (a;) 
\lek{-uj)Al{-uj)Ak{-uj) - -fMiuj)Al{u)Ae{uj)] p'^ = 0, 



ui k/ 

as the sum in uj goes from — Wmax to Wmax for some maximum difference 
between energies Wmax- This proves that is a steady state. Of course this 
does not implies the convergency of any initial state to it (an example is the 
pure dephasing equation derived in section 16.2.51) . Sufficient conditions for 
that were given in the section 15.41 
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6.3 Singular coupling limit 



As explained in section 16. 1[ the singular coupling limit is somehow the com- 
plementary situation to the weak coupling limit, since in this case we get 

— )■ by making the correlation functions to approach a delta function. 
However, as a difference with the weak coupling limit, the singular coupling 
limit is not so useful in practice because as its own name denotes, it requires 
some "singular" situation. For that reason and for sake of comparison, we 
sketch here the ideas about this method from the concepts previously defined. 
For further details the reader will be referred to the literature. 

There are several ways to implement the singular coupling limit, one 
consists in rescaling the total Hamiltonian as p| 1^125] 

H = Ha + a~^HB + 

and perform the "singular limit" by taking a — 0. Starting from the equation 
flS^ under the assumptions of W{t)V = and p(0) = Pa(0) ®Pb(0) we have 

Vp{t) = VpiO) + \ [ ds [ duVV{s)g{s,u)QV{u)Vp{u). (112) 
« io io 

For a moment, let us take a look to the expression at first order in the 
expansion of ^(s, u): 

Vp{t)=Vp{0) + \ [ ds [ duVV{s)V{u)Vp{u) + 0{a-^), 



a' 



which we rewrite as 



PA{t) = Pa{0) - — ds du V{s)[V{u)pa{u)®pb] +0{a-^). 



By introducing the general form of the interaction Eq. ([9T 



PAit) = Pa(0) 



At,{u)pA{u),Akis) 



ds / duCki{s,u) 
Jo L 

+ C*,(s,u) [ifc(s),p^(u)i,(n)] (113) 
here the correlation functions are 



CkAs.u) 



i:i[Bk{s)B,{u)pB\. 



If we assume that p^ is again some eigenstate of the free Hamiltonian we get 

1 



Cm{s - u) 



Tr[5fc(s - u)BepB], 



(114) 
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and by using the eigendecomposition flU^ of and taking into account the 
factor in the free Hamiltonian Hb, we find 



C> 



ke{s 



u 



iuj{s — u) 

due Ti[Bk{uj)BepB]. 



Because the proposition I6.H in the hmit a — ?■ the above integral tends to 
zero as a square integrable function in x = — u). It implies that it 

tends to zero at least as and so we expect this convergence to be faster 
than the prefactor in the correlation functions. However in the peculiar 
case of s = u the correlation functions go to infinity, so indeed in this limit 
they approach a delta function Cki{s — u)(x 6{s — u). Next, by making again 
the change of variable u — )■ s — m, equation flll3p reads 

Pa(^)=Pa(0) + j ds j duCkiiu) Ai{s - u)pa{s - u), Ak{s) 
Jo Jo 

+ C;^{u) [ifc(s), pa{s - u)Aeis - u)] + 0{a-^). 



Therefore in the limit of a — )■ 0, we can substitute m by in the operators 
and extend the integration to infinite without introducing an unbounded 
error 



PA{t) = Pa{0) + I dsTki 
'o 



+ n 



Ae{s)pA{s),Ak{s) 
Ms),pA{s)Ak{s)] 



where 



duCkiiu). 



We may split this matrix in Hermitian j^i = duCke{u) and anti-Hermitian 
part iSke as in equation fllOll) . and differentiate to obtain 



dpAjt) 
dt 



+ 7k£ 



Ae{t)pA{t)Ak{t) - ^ |ifc(t)i,(t),p^(t)} 



where Hhsit) = £^keAk(t)A£(t). Since 'jki has the same expression as 
with = it is always a positive semidefinite matrix and up to 



m 



second order the evolution equation has the form of a Markovian master 
equation. However, from the former it is not clear what happens with the 
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higher order terms; a complete treatment of them can be found in |1U8H110] . 
where it is proven that they vanish. The idea behind this fact is that by 
using the property that for a Gaussian state (which is assumed to be) 
the correlation functions of any order can be written as a sum of products 
of two time correlation functions flll4p . Particularly odd order correlation 
functions vanish because of the condition VV{t)V = 0, and we obtain a 
product of deltas for even order. Actually the expansion is left only with 
those products in which some time arguments appear in "overlapping order" 
such as 

5{ti-t3)5it2-U), 

for ^4 > > ^2 > ^1; and because the time-ordered integration in the 
expansion of Q{s,u) they do not give any contribution. 

Finally by coming back to Schrodinger picture we obtain that in the 
singular coupling limit 



dt 



- i[HA + H^s,pit)] 
+ Ike 



AipA{t)Ak - ^ {AkAe,pA{t)} 



note that i/LS = Ylke ^keA^Ai does not commute in general with H^. 

Despite the singular coupling limit is not very realistic, it is possible to 
find effective equations for the evolution in which self-adjoint operators Ak 
also appear. For example under interactions with classic stochastic external 
fields |108] or in the study of the continuous measurement processes [2T |I125] . 



6.4 Extensions of the weak coupling limit 

Here we will discuss briefly some extensions of the weak coupling limit for 
interacting systems, the reader can consult further studies as for example. 



6.4.1 Weak coupling for multipartite systems 

Consider two equal (for simplicity) quantum systems and their environment, 
such that the total Hamiltonian is given by 

H = Hai + Ha2 + He + aV. 

Here Hai {Ha2) is the free Hamiltonian of the subsystem 1 (2), He is the 
free Hamiltonian of the environment and aV denotes the interaction between 
the subsystems and the environment. This interaction term is assumed to 
have the form V = Xlfc ^i^^ ® -Bfc + Af^ ® Bk, where A\!'^ i^k^^) are operators 
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acting on the subsystem 1 (2), and Bj^ on the environment. So that the 
weak couphng hmit for the whole state of both subsystems pi2 leads to a 
Markovian master equation like 

"^^''^^^ = -^[HAl + HA2 + a'H^s.Pl2m (115) 



dt 



where Ak{uj) are linear combinations of A^^^ and A^^"^ decomposed in eigenop- 
erators of [Hai + Ha2) ■]• Of course the form of these operators depend on the 
form of the coupling. For instance, suppose that we have two identical envi- 
ronments acting locally on each subsystem V = A^j}^ ® B^^ + ■* ® B^^^ , 
then we will obtain 



dt 



UJ k,£ 



A^Pico)p^2m^^^\co) - ^{Al^^\co)A^p{co),p^2{t)} 



Consider now the case in which we perturb the free Hamiltonian of the 
subsystems by an interaction term between them of the form (3Vi2 (where 
P accounts for the strength). Instead of using the interaction picture with 
respect to the free evolution Hai+Ha2+Hb, the general strategy to deal with 
this problem is to consider the interaction picture including the interaction 
between subsystems Hai + Ha2 + PVi2 + Hb, and then proceeds with the 
weak coupling method. 

However if (3 is small, one can follow a different route. First, taking the 
interaction picture with respect to Hai + Ha2 + Hb, the evolution equation 
for the global system is 

j^pit) = -2/3[yi2(t),p(t)] - ia[V{t)rp{t)] = f3Vi2it)pit) + aVimt)- 
Similarly to section 16. H we take projection operators 

Ipp(t) =/5PVi2(t)p(t) + aPV5B(t)p(t), (116) 
^Qp(t) = f3QVi2{t)p{t) + aQVsB{t)m, (117) 
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and find the formal solution to the second equation 

Qp(t) = 6?(t,0)Qp(0) + /3 f dsg{t,s)QVi2{s)Vp{s) 

Jo 

+a [ dsgit,s)QVis)Vpis), (118) 
Jo 

where 

Now the procedure is as follows, we introduce the identity 1 = V + Q just 
in the last term of equation f lll6p . 

Ipp(t) = pVVuit)pit) + aVV{t)Vpit) + aVVit)Qp{t), 
whose formal integration yields 

Vp{t) = Vp{0)+f3 [ dsVVi2{s)p{s) + a [ dsVV{s)Qp{s), (119) 
Jo Jo 

here we have used the condition VV(t)V = 0. Next we insert the formal 
solution (11181) into the last term. By assuming again an initial factorized 
state "subsystems(8)enviroment" (Qp(to) = 0) we find 

Pp(t)=Pp(0) + /3 / dsVVi2{s)p{s) 

Jo 

+ ds dulCi{s,u)Vp{u) + ds dulC2{s,u)Vp{u), 
Jo Jo Jo Jo 

where the kernels are 

/Ci(s,u) = al3VV{s)G{s,u)QVx2{u)V = 
)C2{s,u) = a^VV{s)g{s,u)QV{u)V. 

The first one vanishes because Vi2{s) commutes with V and QV = 0, and 
the second kernel up to second order in a and /3 becomes 

/C2(s, u) = a^VV{s)QV{u)V + C(a^ a^f3) = a^VVis)V{u)V + C(a^ a^P). 

Therefore the integrated equation of motion reads 

Vp{t) = Vp{0) + (5 f dsVVu{s)p{s) 

Jo 

+ [ ds [ duVV{s)V{u)Vp{u) + 0{a\a^^). 
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If we consider small intercoupling /3 such that a > /3, we can neglect the 
higher orders and the weak coupling procedure of this expression (c.f. section 
I6.2p will lead to the usual weak coupling generator with the Hamiltonian part 
/3Vi2 added. This is in Schrodinger picture 

"^^''^^^ = -^[HAl + HA2 + /3V,2 + a'H^s,Pl2m (120) 



dt 



Compare equations (11151) and (11201) . 

Surprisingly, this method may provide good results also for large /3 under 
some conditions, see 



6.4.2 Weak coupling under external driving 

Let us consider now a system A subject to an external time-dependent per- 
turbation PHext{t) and weakly coupled to some environment B, in such a 
way that the total Hamiltonian is 

H = HA + f3H,^,{t) + HB + aV. 

In order to derive a Markovian master equation for this system we must take 
into account of a couple of details. First, since the Hamiltonian is time- 
dependent the generator of the master equation will also be time-dependent, 

whose solution defines a family of propagators £(t2,ti) such that 

Psit2) = £(t2,h)Psiti), 

This family may be contractive (i.e. Markovian inhomogeneous) or just even- 
tually contractive (see section [23j). Secondly, there is an absence of rigorous 
methods to derive at a Markovian master equation (i.e. contractive family) 
in the weak coupling limit when the system Hamiltonian is time-dependent, 
with the exception of adiabatic regimes of external perturbations |126lll27j . 
Fortunately if the ifext(^) is periodic in t, it is possible to obtain Markovian 
master equations, even though the complexity of the problem increases. 
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On one hand, if /3 is small, in the spirit of the previous section we can 
expect that the introduction of PH^xtit) just in the Hamiltonian part of the 
evolution would provide a good approximation to the motion, 

dpA{t) 



dt 



-i[Ha + PHoAt) + a'H^^.pAit)] (121) 



ijj kl 



A,{uj)pA{t)Al{u) - ^{4MA,(a;),p^(t)} 



which has the form of a Markovian master equation (H6l) . 

On the other hand, for larger /3 the only loophole seems to work in the 
interaction picture generated by the unitary propagator 

Taking = without loss of generality, the time evolution equation for 
p[t) = W{t,0)p{t)U{t,0) is 

^ = -za[V{t),p{t)]- (122) 

Following an analogous procedure which was used in section 16.21 for time- 
independent generators, one immediately deals with the problem that it is 
not clear whether there exists something similar to the eigenoperator decom- 
position of V{t) = W{t,0)VU(t,0) as in Note however that since the 
dependency of the operator ifoxt(^) with t is periodical, by differentiation of 
U{t, 0) one obtains a differential problem for each Ak in fl9T]) 



= -z[A,{t),HA + moAt)]- (123) 

with periodic terms. This kind of equations can be studied with the well- 
established Floquet theory (see for example |103[I104] ). Particularly it is 
possible to predict if its solution is a periodic function. In such a case, the 
operator in the new picture will have a formal decomposition similar to (1941) . 
^k(t) = ^fc(w)e*'^*, where now Ak{u) are some operators which do not 
necessarily satisfy a concrete eigenvalue problem. However note that the 
importance of such a decomposition is that the operators Ak{u) are time- 
independent. This allows us to follow a similar procedure to that used for 
time-independent Hamiltonians with the secular approximation, and we will 
obtain 
dpA(t) 



dt 



- i[a i^LS, PA{t)] 

ui kl 



A,iu)pAit)Aliu) - ^{Aiiu)Aeiu),pAm 
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where i^LS = J2ujJ2ke^ki{^)^ki^)^ii^) course. By coming back to 
Schrodinger picture, 

"^P""^^^ -- - z[HA + PH,At) + a'H^s{t),PA{t)] 



dt 



Ae{u,t)pA{t)Al{u,t) - ^{Al{u,t)Ai{u,t),pA{t)} 



here Ak{u},t) = U{t,0)Ak{u)U^{t,0). For these manipulations it is very 
useful to decompose (when possible) the unitary propagator as product of 
non-commuting exponential operators U {t, 0) = ^-I'^ot ^-iHt ^ where the trans- 
formation e*-^o*if(t)e-*-^o* = H removes the explicit time-dependence of the 
Hamiltonian. 

A more systematic approach to this method based on the Floquet theory 
can be found for instance in the work of H. -P. Breuer and F. Petruccione 
[2n[T28]. and in [T29] . 

There is a very important thing to note here. As a difference to the 
usual interaction picture with respect to time-independent Hamiltonians, 
the change of picture given by a time-dependent Hamiltonian can affect the 
"strength" a of ^ in a non-trivial way. If there exist some time, t = 
such that ||ifext(ic)|| = oo, the validity of the weak coupling limit may be 
jeopardized. For example one may obtain things like "1/0" inside of V{t). 



7 Microscopic description: non-Markovian case 

The microscopic description of non-Markovian dynamics is much more in- 
volved than the Markovian one, and developing efficient methods to deal 
with it is actually an active area of research nowadays. One of the reasons 
for this difficulty is that the algebraic properties like contraction semigroups 
and/or evolution families are lost and more complicated structures arise (e.g. 
evolution families which are only eventually contractive). We have already 
mentioned the advantages of being consistent with the UDM description if 
we start from a global product state. According to section 14.31 if the evo- 
lution family £(t,s) describes the evolution from s to t and the initial global 
product state occurs at time to, then it must be eventually contractive from 
that point, i.e. ||i^(t,to)|| < 1 for t > to- However it is possible that > 1 

for arbitrary t and s. Up to date, there is not a clear mathematical char- 
acterization of evolution families which are only eventually contractive and 
this makes it difficult to check whether some concrete model is consistent. 

There are non-Markovian cases where the environment is made of a few 
degrees of freedom (see for example |130H132] , and references therein) . This 
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allows one to make efficient numerical simulations of the whole closed system 
and after that, just tracing out the environmental degrees of freedom, we 
obtain the desired evolution of the open system. 

However when the environment is large, a numerical simulation is com- 
pletely inefficient, unless the number of parameters involved in the evolution 
can be reduced (for example in case of Gaussian states or using numer- 
ical renormalization group procedures |133[I134] ). For remaining situations, 
the choice of a particular technique really depends on the context. For con- 
creteness, we sketch just a few of them in the following sections. 

7.1 Integro-difFerential models 

As we have seen in section 16.11 the exact evolution of a reduced system can 
be formally written as an integro-differential equation ( 186|) . 



Typically this kernel cannot be written in a closed form and, furthermore, 
the integro-differential equation is not easily solvable. On the other hand, it 
is possible to perform a perturbative expansion of the kernel on the strength 
of V{t) and proceed further than in the Markovian case, where the series 
is shortened at the ffist non-trivial order (for an example of this see [135] ). 
Of course, by making this the complete positivity of the reduced dynamics 
is not usually conserved; and one expects a similar accuracy for the time- 
convolutionless method (at the same perturbative order) explained in the 
next section, which is simpler to solve. 

As an alternative to the exact integro-differential equation, several phe- 
nomenological approaches which reduce to Markovian evolution in some lim- 
its has been proposed. For example, if C is the generator of a Markovian 
semigroup, an integro-differential equation can be formulated as 



where k{t — t') is some function which accounts for "memory" effects, and 
simplifies the complicated Nakajima-Zwanzig kernel. Of course, in the limit 
k{t—t') — )■ 6(t—t') the Markovian evolution is recovered. The possible choices 




with a generally very complicated kernel 



}C{t,u) =VV{t)g{t,u)QV{u). 




(124) 
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for this kernel which assure that the solution, £(t,o), is a UDM have been 
studied in several works |136H138] . being the exponential ansazt k{t — t') oc 
ge^i^-^') the most popular. 

Another interesting phenomenological integro-differential model is the so- 
called post-Markovian master equation, proposed by A. Shabani and D. A. 
Lidar 

^ = C J\it - t')e^^'-'\it')]. (125) 

It also approaches the Markovian master equation when k{t — t') — )■ S{t — t'). 
Further studies about the application of these models and the conditions to 
get UDMs can be found in [TiOHn2] . 

Apart from these references, the structure of kernels which preserve com- 
plete positivity has been analyzed in [143] . 

We should note however that, it has been recently questioned whether 
these phenomenological integro-differential equations reproduce some typical 
features of non-Markovian dynamics |144j . 

7.2 Time-convolutionless forms 

The convolution in integro-differential models is usually an undesirable char- 
acteristic because, even when the equation can be formulated, the methods 
to find its solution are complicated. The idea of the time-convolutionless 
(TCL) forms consist of removing this convolution from the evolution equa- 
tion to end up with an ordinary differential equation. Behind this technique 
lies the property already expressed in section H75] that a UDM, £{t,to)i can be 
expressed as the solution of the differential equation (l40l) 

dp Ait) 



with generator Ct = — ^jir^£(t\ v general this equation generates an evolu- 



dt ^(t,to)' 

tion family that is not contractive, i.e. ||i^(t2,ti) || > 1, for some ti and t2 and 
so non-Markovian. However, it generates an eventually contractive family 
from to, ||^^(t,to)ll — 1; starting from to, £{t,to)^ usually assumed 

to be a UDM by hypothesis. 

To formulate this class of equations from a microscopic model is not 
simple, although it can been done exactly in some cases such as, for a damped 
harmonic oscillator |145yi46j or for spontaneous emission of a two level atom 
|147pi48] . In general, the recipe is resorting to a perturbative approach again. 
Details of this treatment are carefully explained in the original work |149j 
or in some textbooks (see [201 [21] instance). The key ingredient is to 
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insert the backward unitary operator V({s,t) where s < t defined in 
p(s) = U{s,t)p{t), into the formal solution of the Qpit), Eq. 



Qp(t) = ^?(t,0)Qp(0)+ / dsg{t,s)QV{s)VV({s,t)p{t). 

Jo 

Now, by introducing the identity 1 = V + Q between U{s,t) and p(t), and 
after considering a factorized initial condition we obtain 

[l-T(t)]Qp(t) = T(t)Pp(t), 

here ^ 

T{t) = a [ dsg{t,s)QV{s)VU{s,t), 



where we have rewritten V as aV. 

Next, note that T(0) = and T{t)\a=Q = 0, so that the operator [1 — T(t)] 
may be inverted at least for not too large couplings and/or too large t. In 
such a case we would obtain 

Qpit) = [l-T{t)]''T{t)V pit). 



Introducing this in the second member of equation (1791) . the first term van- 
ishes as W{t)V = 0, and we arrive at the differential equation 

dVpjt) ^ 

where 

= aVVit)[t - T(t)]-iT(t)P. 
We can formally write the inverse operator as a geometric series 

oo 

[l-Tit)]-' = J2T-it), 

n=0 

and then the generator as 



= aVVit) J2 



n=l 



On the other hand, T"(t) can be expanded in powers of a by introducing the 
expansions of ^(t,s) and Uis,t), thus we can construct successive approxi- 
mations to Ct in powers of a, 



;(n) 

n=2 
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Here the sum starts from n = 2 because at lowest order T(t) ~ a dsQV{s)V, 
and 

nt nt 

£f ) = / dsVV{t)QV{s)V = / dsVV{t)V{s)V, 
Jo Jo 

which leads actually to the differentiated version of equation (190|) . as expected 

at lowest order. After a little bit of algebra, the next term in the expansion 

turns out to be 

£f) ^ dtiVV{t)V{ti)V{s)V, 

and further terms can be calculated in the same fashion. 

In addition, D. Chruscihski and A. Kossakowski jl50j have recently shown 
the equivalence between integro-differential equations and TCL equations. 
They connect both generators by means of a Laplace transformation showing 
some kind of complementarity between both methods. That is, if the integral 
kernel is simple, the TCL generator is highly singular, and viceversa. 

On the other hand, given a TCL equation and an initial condition at 
time to, the problem of knowing whether its solution S(^t,to) is a UDM, i.e. a 
completely positive map, is an open question for the non-Markovian regime. 
Note that the perturbative expansion at finite order of the TCL generator 
can generally break the complete positivity (in the same fashion as a finite 
perturbative expansion of a Hamiltonian dynamics breaks the unitary char- 
acter). So the classification of TCL generators which are just eventually 
contractive from a given time to is an important open problem in the field 
of dynamics of open quantum systems. Partial solutions have been given in 
the qubit case [151] and in the case of commutative generators £f |1521I153] . 
For this latter one [£t^,£t2] = 0, for all ti and t2, so that the time-ordering 
operator has no effect and the solution is given by 

Therefore, certainly £(t,to) is completely positive if and only if J^* Cfdt' has 
the standard form (1511) for every final t. 

7.3 Dynamical coarse graining method 

The dynamical coarse-graining method was recently proposed by G. Schaller 
and T. Brandes |154[ll55j . Thiswas suggested as away to avoid the secular 
approximation under weak coupling while keeping the complete positivity in 
the dynamics. Later on, this method has been used in some other works 
|156[I157] . Here we briefly explain the main idea and results. 
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Consider the total Hamiltonian f lTT]) and the usual product initial state 
p(0) = Pa(0) ® Pbi where Pb is a stationary state of the environment. In the 
interaction picture the reduced state at time t is 



PA{t) = Tie [U{t, O)p^(O) ® PbU^I, 



(126) 



where 



is the unitary propagator, and a again accounts for the strength of the in- 
teraction. Next we introduce the expansion of the propagator up to the first 
non-trivial order in equation ^1 



pA{t)= Pa{0)-y'^ I dti I dt2TTB 



V{ti), V{t2),PA{0)(^PB +0{a' 



here we have assumed again that the first order vanishes Tr[y(t)p^] = like 
in (185|) . Let us take a closer look to structure of the double time-ordered 
integral 



r dU dt2TTB V{h), V{t2),PA{0)®PB 



dti / dt2d{h - t2) TiB \v{h), \v{t2),PA{^) ® PB 



+ 



dh / dt2e{t2-ti)TlB V{t2), V{ti),pA{0)^pB 



cf\pAm+^f[pAm+^T[p 



,(2) 



(2), 



where by expanding the double commutators we have found three kind of 
terms, £p\ £3^^ and £3^^ The first one appears twice, it is 



(2) 



-[Pa(0)] 



- dh «(ti-t2)TrB Vih)pA{0)®PBV{t2: 
Jo Jo '- 

+9{t2 - h) Tib [v^(ti)PA(0) ® PBV{t2) 

-[ dh [ dt2[e{h - 12) + e{t2 - h)] 

Jo Jo 

Tib \v{h)pA{0)®pBV{t2] 



- dt, dt2TiB Viti)pAiO)(S)PBVit, 
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because obviously the step functions cancel each other. Taking in account 
the other analogous term (that is the corresponding to interchange ti -H- ^2 
in the double commutators) we can write the total contribution as 



-2TrB[ApA(0)®pBA] 



where A = /J V{t')dt'. 

The second term goes like 

4'^Ipa(0)] = [ dh [ dt2e{h-t2)TTB\v{h)V{t2)pA{0)(^PB 
Jo Jo 

+9{t2 - h) Tib [v{t2)V{h)pA{0) ® pb 

= [ dti [ dt2[e{h - ta) + 0{t2 - h)] 
Jo Jo 

Ttb [v{t2)V{h)pA{0)^pB 

+9{t^ - t2) Trs {[V^(ti), V{t2)]pA{^) ® Pb 
= Ttb[A^Pa{0)®Pb] 

+ ^ dt2eih - t2)TrB {[V{h),V{t2)]pAiO) ^ Pb] . 

And similarly the remaining term can be expressed as 
^?[PAm = Tr^ [pa(0)®PbA2] 

- j^dhj^ dhOih - t2) Trs {pa(0) ® PB[V{ti), V{t2)] 
Thus, everything together gives pAit) = Pyi(O) + £*[pyi(0)] + 0{q^) with 



a 



^ [o^\PAm + oi\PAm + cf{PAm 



:f)[PA(o)] + £r^ 

^ -2[i7*s,PA(0)]+I?*[pA(0)]. 
Here the self-adjoint operator f^Ls 



2 rt 



His = Yi "^^^ '^^^^(ti - t2) Tie \ [V{h),V{t2)]pB \ , 
and the "dissipative" part is given by 
P*[pA(0)] = «2TrB 



1 



Apa(O) ® pbA - - (A^, p^(0) ® Pb} 



88 



By using the spectral decomposition of it is immediate to check that for 
any fixed t, £* has the Kossakowski-Lindblad form flM]) . 
Since at first order 

PA{t) = (1 + £*)pa(0) + 0{a^) - e^VA(O), (127) 

the "dynamical coarse graining method" consists of introducing a "coarse 
graining time" r and defining a family of generators of semigroups by 




And the corresponding family of differential equations 

p^(t)=£>~^(t), 

which solution 

p\{t) = e*^">A(0), (128) 
will obviously coincide with f ll28p up to second order when taking t = t, this 

is p-^=\t) = pA{t). 

Therefore we have found a family of semigroups fll27p which have the 
Kossakowski-Lindblad form for every r and give us the solution to the evolu- 
tion for a particular value of r (r = t). This implies that complete positivity 
is preserved as we desire. In addition, one may want to write this evolution 
as a TCL equation, i.e. to find the generator Ct such that the solution of 

j^PAit) = CtiPAit)] 

fulfills pA{t) = p\{t). That is given as in fjlO|l . 



Ct = 



—e'- 


e ^ = 


d rt 


_dt 




_dt 



Particularly for large times t — )■ oo, this generator Ct approach the usual 
weak coupling generator (11041) (in the interaction picture). This is expected 
because the upper limit of the integrals in the perturbative approach will 
approach to infinity and the non-secular terms e^^^ will disappear for 
large t (because of proposition 16. ip . This is actually equivalent to what we 
obtain by taking the limit a — )■ in the rescaled time. Details are found 
in [T5i] . 
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8 Conclusion 



We have reviewed in some detail the mathematical structure and the main 
physical features of the dynamics of open quantum systems, and presented 
the general properties of the dynamical maps underlying open system dynam- 
ics. In addition, we have analyzed some mathematical properties of quantum 
Markovian processes, such as their differential form or their steady states. 
After that we have focused our attention on the microscopic derivations of 
dynamical equations. First of all we have analyzed the Markovian case, re- 
calling the weak and singular coupling procedures, as well as some of their 
properties; and secondly, we have discussed briefly some methods which go 
beyond Markovian evolutions. 

Several interesting topics have not been tackled in this article. Examples 
are influence functional techniques |21[|23[l33yi58yi59] . stochastic Schrodinger 
equations [T5 | [26 t[T60l - ll63j . or important results such as the quantum regres- 
sion theorem, which deals with the dynamics of multi-time time correlation 
functions [T ^ 12 ^ 1164^1165] . Moreover, recently several groups have ad- 
dressed the important problem of introducing a measure to quantify the 
non-Markovian character of quantum evolutions |166H173] . However, despite 
the introductory nature of this work, we hope that it will be helpful to clarify 
certain concepts and derivation procedures of dynamical equations in open 
quantum systems, providing a useful basis to study further results as well as 
helping to unify the concepts and methods employed by different communi- 
ties. 
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